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We observe that the large world sheet RG in c = 1 matrix model, formulated in 0|2], with A^^ 
quantum mechanical degrees of freedom at small compactification radius is capable of capturing 
dimensional mutation. This manifests in deforming the familiar AdS2 quantum mechanics in 
the minisuperspace Wheeler-de Witt (WdW) cosmology of the 2D quantum gravity, obtained 
by the large N RG with quantum mechanical degrees of freedom only, to a modified WdW 
cosmology describing tunneling to an inflationary de Sitter vacuum and its evolution. The scale 
fluctuation plays an important role in providing an ansatz for uniquely choosing the initial 
wave function. We observe that the nonperturbative effects due to the A^^ quantum mechanical 
degrees of freedom introduce explicit open string moduli dependence in the wave function via 
the Hubble scale, which determines the geometry of the true vacuum one tunnels to. The 
modified WdW equation also captures controlled formation of baby universes of vanishing size 
that self-tunes the de Sitter cosmological constant to be small positive. 
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1 Introduction 

One of the central problems in quantum grayity in recent time is to find a yacuum waye function 
for a consistent quantum mechanical description of the early uniyerse (with respect to the 
recent obseryational data) and to choose it uniquely, if there is a dynamical selection principle. 



2 



A simple quantum cosmological approach [3] (see j3] for a review) would be to arrive at a 
Schrodinger like (WdW) equation that would describe tunneling of the universe from nothing 
to an inflationary de Sitter vacuum rather than to an Euclidean Hartle-Hawking vacuum 
[H]. Once the schrodinger equation is derived, one would also like to have its unique solution 
[7j. We refer to [1111111011111112111311111113 for a discussion on quantum creation of inflationary 
universe from nothing. One can see recent works, for example JH] and references therein, for a 
comparative discussion of the two types of wave functions in string theory context. 

In string theory the universe can be chosen to be the holes in the world sheet, that can be 
filled with D-branes. The dimensionality of the universe and the issue of openness/closedness 
or the homogeneity /inhomogeneity depends on how we fill these holes. Even though these holes 
are world sheet quantities, they can have a good space-time description. The one point function 
of these objects describes the wave function ip of the universe with a scale factor / = j e^'^'^(Pa 
[if is the Liouville field in a noncritical theory). In minisuperspace approximation, the WdW 
equation (a BRST constraint of the field theory) 

i(Lo + Zo)|^) = (1.1) 

is solved for Euclidean Hartle-Hawking vacuum with homogenous FRW cosmology. Because of 
reaching a huge size in a small time interval, for most of the purposes the cosmological evolu- 
tion of the inflationary universe is semiclassical and minisuperspace computations successfully 
capture a lot of features PHj- However, in an inflationary universe, the long wave fluctuations 
of the scalar field drives geometry at much larger scales (at / > H^^ ~ M^^) to be highly inho- 
mogeneous, far away from that of a homogeneous Friedmann space. The quantum fluctuation 
of the scalar field may drive some parts of the universe to inflate eternally, whereas other parts 
pass the end of inflation and undergo different phases of evolution. Much of the important 
inflationary physics are actually lost in the realm of the minisuperspace approximation, that 
can not get rid of an initial presumption of global homogeneity of the universe. The initial 
wave function in early universe dealing a cosmological singularity with particle production on 
the horizon is essentially quantum mechanical and is worth seeking a proper stringy description 
(see 120] for a review on the subject; also see [21] for a recent study of decay of metastable vacua 
in 2D Liouville gravity, for issues on nucleation in 2D Liouville gravity). 

The problem can be defined in the following way. The cosmological singularities often 
can be resolved by a phase of suitably chosen closed string tachyon condensate much like the 
exponential wall in a time-like Liouville theory (see [221123] and references therein). Generically 
the fields in such tachyon phase become heavy and tend to evacuate configurations of real or 
virtual particles or excitations that source any component of string field (see [23] and further 
references). More specifically the S"- matrix of the system, that governs the perturbative BRST 
consistency condition on the states in the tachyon condensate phase, does not have perturbative 
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asymptotic particle poles. As a result string world sheet reaches spatial infinity in a finite world 
sheet time, giving rise to holes in the world sheet. Though in special cases D-brane boundary 
states make these holes consistent, generically they give rise to BRST anomaly that can be 
resolved by considering correlations of these holes corresponding to their unitary mapping. 
These correlations can be nonlocal representing interaction between multiple world sheets [23 
|2ni and can be complicated enough to model the dynamics in the black hole interior (see 
examples discussed in [21] )• Note that the above scenario of preserving unitarity is reminiscent 
of the black hole final state proposal for solving the information paradox j27j. However, imposing 
cancelation of BRST anomaly alone does not specify a unique correlation and thus a unique 
final state, rather yields a restricted set of states in the tachyon condensate phase. 

In the light of above discussion the vacuum wave function for the early universe can be 
described by a wave function for some kind of tachyon condensate (T) that replaces the initial 
cosmological singularity. A simplest undeformed condensate would correspond to the Euclidean 
Hartle-Hawking vacuum whose evolution is given by the usual WdW equation ^. However, 
the stringy modification to the WdW equation can come from an anomaly in the BRST like 
constraint that is resolved by a nontrivially deformed condensate (T + 5T) (See jHIIl IHH 
122 E21 for other relevant approaches/proposals of modifying the Hartle-Hawking vacua). The 
allowed set of states surviving in this nontrivial deformed tachyon condensate phase can arise 
from various 'beyond minisuperspace' contributions. As we discussed before, they can arise from 
nonlocal correlations of the holes via multi-trace [23 12^ or topology changing amplitudes that 
are essentially beyond minisuperspace [3]. These amplitudes describe emission or absorption 
of baby universes in the WdW cosmology [33J [SHI (HHl HZ] that can self-tune the cosmological 
constant to appropriate value. However, in a full nonperturbative setup (like the one we are 
going to use in this paper) the deformed condensate can contain nonperturbative objects that can 
feel all the nonzero modes of dilaton. Note that, in the perturbative computation of Euclidean 
vacuum in [2H|, the tachyon condensation effectively masses up all the closed string modes 
and hence the fluctuations in dilaton too, freezing it to its bulk weakly coupled value. Thus 
the setup remains perturbative there. There are other similar evidences of perturbativity in 
the closed string tachyon condensation such as repulsion felt by D-brane probe from winding 
tachyon phase [HE] (also see fate of twisted D-branes in orbifolds liH]). 

However, in this work we will argue that the presence of nonperturbative objects, that 
penetrate the tachyon wall down to strong coupling (also see the 2D example in [33), make 
it possible to feel the dilaton fluctuations which play a vital role in uniquely determining the 

^In time dependent background in perturbative string theory (in type IIB and as well as in heterotic strings) 
such a condensate is constructed from the correlators supported in the weakly coupled bulk j^H] that seems to 
be a useful setup to study space-like big bang singularities in Milne space-time. For a nonperturbative study, 
see [2ni- 
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quantum state corresponding the nontrivially deformed tachyon condensate. We will also discuss 
how the ability to feel the dilaton fluctuations can have implications in resolving a cosmological 
singularity. In particular we note the possibility of using the properties of density perturbation 
spectra, such as its near flatness, to determine an ansatz for dynamical vacuum selection via 
the initial profile of the scale (dilaton) fluctuation. 

In the language of WdW cosmology, computing the wave function for the appropriate de- 
formed tachyon condensate replacing a cosmological singularity would require going beyond the 
minisuperspace by considering general matter dependent deformations and deformations due 
to higher order derivatives in scale. In other words one needs to access the full superspace by 
considering the arbitrary scale and matter fluctuations due to all the nonzero modes of dila- 
ton and matter. The nonzero modes of matter can be excited by appropriately considering 
Neumann boundary condition on the world sheet holes. However, the role of arbitrary matter 
fluctuations in a nonperturbative tachyon condensate we are considering, is an interesting issue 
that we take up elsewhere Thus studying an initial wave function for an inflationary vac- 
uum eventually would be tied to a time dependent nonperturbative setup capturing dynamical 
fluctuations and nonlocal processes. Even though we start with a (sub) critical background, we 
believe the phenomena we are observing is generic due to dimensional mutation brought by 
these nonperturbative objects, that we explain below. 

So far we focussed on the question of capturing the inflationary nature of the initial wave 
function and its unique solution. Now the issue of getting a (Lorentzian) de Sitter vacuum as 
opposed to an (Euclidean) AdS is tied to the question of accessing the supercritical regime of 
the theory instead of a (sub)critical one [33 EHl EZI- Let us take up the example of bosonic 
string theory. In 2D noncritical bosonic string theory (Liouville coupled with c = 1 matter) 
the minisuperspace WdW equation for (T), both in the continuum and discrete version (free 
fermionic picture of matrix quantum mechanics) is a Bessel equation solving for an Euclidean 
AdS2 vacuum {ip ~ Kq^^/JII) , ^ > 0) (see the reviews and [IHI)- The minisuperspace 

equation in a supercritical setup (Liouville coupled to c > 25) with a tachyon condensate 
in one of the matter directions solves for global de Sitter in the far future {Hq{^—iiI)) ^6] . 
However, due to the simplicity of the condensate, the vacuum wave function does not capture 
any kind of inflationary scenario or proliferation of the de Sitter. The striking feature is that, 
the minisuperspace WdW equation in this case is apparently a de Sitter quantum mechanics of 
single matter but with other c = 25 matter directions fixing the Weyl anomaly. Thus their only 
role is in fixing the de Sitter sign of the cosmological constant, the Lorentzian signature of the 
world sheet metric and a sensible interdependence of the scale and the string coupling for a viable 
cosmological scenario. On the other hand, remnants of dS2 quantum mechanics is already there 
even in minisuperspace 2D noncritical theory. This is evident from the free fermion quantum 
mechanics, where a special class of time (or c = 1 matter) dependent Wqo deformations of the 
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Fermi sea gives rise to past and future space-like boundaries like a space-time cosmological event 
horizon |17] ^. Let us here remind ourselves that the semi classical condensate (T) giving rise to 
Euclidean Hartle Hawking vacuum {AdS2) of the c = 1 theory, can be thought of as simplest 
fluctuations on the Fermi surface corresponding to an inverted harmonic oscillator potential 
|43j . Thus perhaps in the nonperturbative setup, the nontrivial time dependent deformations 
of this Euclidean Hartle Hawking vacuum in the (sub)critical theory is carrying out some kind 
of dimensional mutation by capturing the physics of a supercritical regime (see discussions 
in jSJ and references therein for comments on dimensional mutation in closed string tachyon 
condensation in the context of cosmological singularities). 

In this paper we put together all the ingredients of the problem in the context of the large 
N RG E], a Wilsonian world sheet RG, of Dirichlet boundaries in 2D noncritical string 
theory (the FZZT branes [H2] with Dirichlet boundary condition on matter). They describe 
the quantum cosmology of the one dimensional universes with homogeneous matter. The wave 
function of the universe is given by the one point function of the loop operator W{l,t) = 
-i:Tr exp[/$(t)] [531 EH mi iSnj- Note that the loop operators we study also have space-time 
interpretation. Here we compute the wave function by performing the large RG of the one 
point function of the loop operator of the noncritical theory around a c = 1 fixed point at a 
string scale compactification radius (a self-dual radius R ~ 1). In fact, in our RG scheme, the 
Dirichlet boundary condition on matter at a c = 1 fixed point prefers a small compactification 
radius. As a result the large RG deals with all the A^^ quantum mechanical degrees of 
freedom, instead of the A^ eigenvalues only as in minisuperspace free fermionic description. 
This translates into unavoidable large fluctuation of dilaton and nonlocal effects associated to 
the presence of the off diagonal elements (that can not be integrated out in c = 1 theory on 
a small circle) and presumably release of large number of nonperturbative states |SZ1 EHl (Hnj 
causing a dimensional mutation. Thus the RG captures a nontrivially deformed condensate that 
modifies the WdW equation by scale derivatives of all orders, a nontrivial time dependence and 
nonlocal deformations corresponding to emission and absorption of point-like baby universes. 
Strikingly the process of emission and absorption of babies are controlled by 0(1/A^). The short 
distance wave function of this modified WdW gives a tunneling wave function that corresponds 
to birth or nucleation of the universe into an inflationary de Sitter vacuum. Here let us briefly 
mention that in the recent semiclassical study of the decay of weakly metastable vacua in 2D 
Liouville gravity in [25 , gravitational inflation inside the region of lower energy phase due to 
droplet fluctuation in the critical swelling process modifies the standard exponential suppression 
of the nucleation probability to a power law one, much like the behavior of the dynamical lattice 
Ising model. Now the far future wave function gives a global de Sitter as would have been given 

^See 0]0HlEni for the issue of particle production and thermality. Also see (SO) for other time dependent 
solutions of similar Woc deformed condensates with a cosmological interpretation. 
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by a dS minisuperspace theory. The scale derivatives are playing an important role. Truncating 
them to quadratic power gives only the global de Sitter (Jo(V~/^0) ^^e far past like a de Sitter 
minisuperspace theory instead of a tunneling wave function. The arbitrary scale fluctuations 
help to fix the initial condition uniquely. 

It is interesting to observe that, the time dependent nonperturbative deformations in non- 
critical theory is doing something nontrivial to enhance the effective space-time dimensions and 
hence accesses the Lorentzian de Sitter geometry dwelling in a supercritical regime. Presum- 
ably the dimensional mutation is due to large number of nonperturbative states belonging to 
the nonsinglet sector of the matrix quantum mechanics appearing at the self-dual radius. The 
real nature of these states (whether they are best described as world sheet vortices or long 
winding strings) is yet to be understood jSZl EHl 1^01 inDj (also see j21 EI] for a discussion on 
the nature of these states in the context of 2D black hole). It has been pointed out in [ST] 
that in the case of spacelike singularities studied in |2H] and in similar cases with positively 
curved spatial slices, the singularity is lifted by a tachyon phase with fewer degrees of freedom, 
as opposed to the cases corresponding to negatively curved spaces. This is a candidate big bang 
case where with forward time evolution along the larger direction of the compactification circle 
the dimensionality (cg//) increases. This is similar to the realistic inflationary cosmology where 
universe evolves in the direction of the evolution of de Sitter entropy [HSinS])- Note that, in our 
case the universe evolution is not in the direction of decreasing R, instead is in the direction of 
increasing /. The big bang singularity in our picture sits at small scale (/ ~ H~^) and is resolved 
by popping up of a simple tachyon condensate (Euclidean vacuum) at large R. Whereas tuning 
the radius to string scale deforms the condensate setting in dimensional mutation to supercrit- 
ical regime. In this case an inflationary de Sitter vacuum pops up at / ~ H^^ to resolve the 
singularity. Nevertheless the universe evolution in our picture does take place in the direction of 
evolving de Sitter entropy. This is because for / < we have tunneling to the corresponding 
vacuum from nothing. Thus in effect the R direction maps to the set of initial condition in 
string theory that maps to the available set of vacuum {AdS or dS) through the large N RG 
(see the figure H below). However, as we will see, the dynamical vacuum selection here depends 
on other factors (dilaton fluctuation) that uniquely specify the initial tunneling wave function. 
As pointed out in [01], the universe we live in seems to have some mechanism of self tuning. 
Perhaps the large RG incorporating dilaton fluctuation and nonlocal effects due to multiple 
world sheets inherently exhibits such self tuning that enables one to flow to a more realistic 
vacuum, starting around a c = 1 fixed point. It will be interesting to generalize the problem, 
in particular for an inhomogeneous universe with arbitrary scale and matter fluctuation, which 
we will report later [32] • 

The plan of the paper is as follows. In section 2 we review and elaborate the problem of 
inflationary de Sitter in 2D cosmology and summarize our results and critical observations. 
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Figure 1: The mapping of the set of initial conditions in string theory to the set of available 
vacua. The broken arrow denotes possible time dependent deformations of the Fermi surface 
picture that may connect the two sets of vacua. 

In section 3 we review our previous work ^ on large RG analysis in c = 1 matrix model. 
In section 4 we deduce the familiar minisuperspace WdW equation of free fermions, namely 
the AdS2 quantum mechanics, from large A^ RG. In section 5 we compute its deformation, 
the generalized WdW equation for one dimensional homogenous universe with arbitrary shape 
and scale fiuctuations. In section 6 we identify the Hubble scale and then solve for the wave 
function at small and large scale limits and discuss the cosmological implications: tunneling to 
infiationary de Sitter in the far past and an expanding global de Sitter in the far future. In 
section 7 we discuss the topology changing amplitudes captured by large A^ RG and suppression 
of baby universes and self-tuning of cosmological constant. 



2 2D cosmology and inflationary de Sitter vacuum 

Let us now review and elaborate on the basic philosophy of the problem in bosonic string theory. 

2.1 The minisuperspace WdW cosmology in (sub) critical and super- 
critical theory 

Let us consider Euclidean quantum gravity described by Liouville field theory coupled to c = 1 
matter or its discrete matrix quantum mechanics version, where one has a good nonperturbative 
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description of the time dependent physics. The tachyon field T{ip, X{a)) on the world sheet is 
directly related to the macroscopic loop fields W{l,t) EH 1^ 1^ 1^ 1^ via 



~ e-5«'^ T((^,t). (2.1) 

Here, if is Liouville field and Q is the slope of the linear dilaton background (Q = \p2 for c = 1 
theory). The relation can be intuitively understood by considering the first order fluctuation 
that corresponds exactly to the tachyon wave function and satisfies minisuperspace WdW 
equation up to a factor of coupling constant. The vacuum expectation (iy(Z,t)) is given by the 
Hartle-Hawking type wave function ^/'(Z). In the matrix quantum mechanics, these are nothing 
but the FZZT branes or the one dimensional boundaries of length / 

V^(/) = l(Tre'*W), (2.2) 
with Dirichlet boundary condition on the world sheet matter 

W{1, tUat 5{j e^'^ - /) 5{X{a) - t) . (2.3) 

Here W{l,t)cont represents the loop operator in the corresponding continuum Liouville theory. 
Their wave function, computed by the eigenvalue representation or the free fermion representa- 
tion of matrix quantum mechanics [331123 EH]; satisfies the same minisuperspace WdW equation 
as in the continuum Liouville field theory (accompanied by proper identification of the matrix 
model couplings with the continuum parameters) 

- (/oj-)' + + ^^h{h)h=o = , /o = e^^« . (2.4) 

OLq J 

In matrix model computation, the all genus version of ()2.4|) for zero Fourier mode of matter 
modifies the WdW by adding a —qIIq term in the potential. Note that, the free fermion 
computation we mentioned here is computed in minisuperspace {i.e. ip space is truncated to 
Lpo). This is manifest in the massless scalar (the fiuctuation of the eigenvalue density) quadratic 
action that exactly maps to the minisuperspace WdW action of the Liouville field in a suitable 
coordinate in the space of eigenvalues j3H] • Both the continuum and the matrix model solutions 
for the minisuperspace Hartle-Hawking vacuum (an undeformed tachyon condensate) are just 
Bessel functions with simple FRW cosmology. The solution decaying at large length is the non- 
normalizable wave function Ky{2^—jjilQ). In fact Kq{2^—iiIq) describes AdS vacuum (/i > 0) 
that has simple big bang and big crunch [3^1 • Thus the free fermion quantum mechanics ()2.4|) 
with /i > is a quantum mechanics for AdS2- From the point of view of loop operator, it 
describes dynamics of spherical loops with a homogeneous matter. Note that a wave function 
given by Bessel Jo(2^— yuZ) and Hankel Hq{2^—^1) would have described the behavior of global 
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de Sitter in the far past and far future respectively |l6], as it will arise in a different context. 
However, as we will see, these Bessel type of solutions describe merely homogeneous FRW 
geometries and do not seem to capture any inhomogeneity introduced by inflationary scenario 
with cosmological time evolution. 

In terms of the free fermion phase space, the undeformed tachyon condensate can be de- 
scribed by the physics of a static Fermi sea corresponding to the eigenvalue hamiltonian with 
inverted oscillator potential filled up to a level fi. The tachyon wave function can be visualized as 
the fluctuations 6W along the Fermi surface H{p, A) = fi. Now a general potential V{X) would 
spatially deform (fold, shift or open) the Fermi surface H{p, A) = — = /i to a Fermi surface 
of an arbitrary shape H{p, A) = p"^ + ^V{X) = fi. This introduces higher derivatives in the 
minisuperspace WdW equation that modifies the genus zero tachyon wave function according 
to 

+ llV{d/dlo) + ^jv'{d/dlo)] 6Wh=o = 2fill 6Wh=o • (2.5) 

Explicit time (matter) dependent deformation (tearing into droplet, opening and draining to the 
other side) of the Fermi surface result in nontrivial time dependent coefficients in the kinetic and 
potential terms of the WdW equation. Effectively it boils down to a time dependent fi{t) and a 
time dependent over all multiplication factor in the tachyon wave function. A class of integrable 
(woo) time dependent deformation have been shown to have cosmological interpretations (e.g. 
bounce cosmology) jHO]. 

For inflationary scenario, one needs an appropriate temporal and spatial deformation of 
the condensate such that the modified WdW gives rise to some kind of nucleation to a de 
Sitter vacuum at UV (or in the far past) and its proliferation. Note that from ()2.4|) . a naive 
continuation to the de Sitter geometry (// < 0) by dialing the sign of cosmological constant fi 
is subtle and apparently does not make sense for the eigenvalue dynamics of the c = 1 matter. 
There the /i — > corresponds to instability due to spilling of the eigenvalues from the top of the 
inverted oscillator potential (the continuum limit in the double scaling procedure) and the fj, < 
side is unbounded. However, considering an interesting class of u^oo deformations corresponding 
to special class of time dependent perturbations of the Fermi sea (opening and draining of the 
Fermi sea to the other side), one can show formation of spacelike horizons either or both at 
past and future spacelike infinities [Hj and also particle production |1HI EH EHj- Strikingly 
these properties indicate remnant of some kind of de Sitter quantum mechanics. Thus it is not 
meaningless in c = 1 matrix model to look for most general (integrable and nonintegrable ^) time 
dependent deformation of the tachyon condensate that will give rise to wave function exhibiting 
nucleation to de Sitter vacuum and possibly its proliferation too. For c > 25 matter a global 



■^AU the integrable Woo deformations can be written in minisuperspace variable (po whereas the most general 
deformations in (accessible) superspace are all not necessarily integrable. 
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de Sitter geometry at large scale can be probed by considering a Liouville-Sine-Gordon model 
of a inhomogeneous tachyon condensate jlHl- The condensate is in the form of gravitationally 
dressed cosine matter potential in one of the matter (spatial) direction and is homogenous in 
remaining scalar fields that only play a role through the Weyl anomaly (they are all massless 
and remain in their initial state throughout the evolution of the universe) . The minisuperspace 
WdW equation is then modified by a matter dependent cosine term in the potential. This gives 
a solution exhibiting global de Sitter geometry at large scale though it recoUapses at small scale. 
However, such a condensate is still too simple to catch proliferation of the de Sitter or particle 
production. 

In this respect one has to keep in mind that the quantum cosmology for c = 1 matter has few 
subtleties that are better dealt in c > 25 regime jHEl EZl EH] • This is because the length scale of 
the universe behaves as / = ^ d?a e'*''^ , where d?ae"''^ is conformally invariant {i.e. (5 = 7 + 2/7). 
For c = 1, Q and hence 7 are real and the length scale grows proportionally with the string 
coupling Qs ~ e^*^*^. The universe is then small at weak coupling and large at strong coupling, 
which is unrealistic for doing cosmology. For c > 25, Q and 7 are purely imaginary. Thus one 
Wick rotates ip — * —iip getting a length scale inversely varying with string coupling Qs ~ e~2*^'^. 
The wick rotation has an added advantage that it makes the kinetic energy of the time-like 
Liouville coordinate negative (analogous to that in Ad Euclidean gravity) so that the world 
sheet geometry has a Lorentzian interpretation. Thus for technical reasons we discussed so far, 
c > 25 matter has been preferred over the c = 1 matter theory in probing de Sitter cosmology 
on the world-sheet jlHl- Finally for c = 1 theory at strong coupling, formation of baby universes 
of macroscopic size is favored that may absorb the inflating core of the universe and force it 
to recoUapse, once the fluctuation drives the infiaton off its potential maxima. In single closed 
universe cosmology of gravity coupled to c = 1 matter in 7?. x S*^, where the winding strings in a 
Sine-Gordon Liouville setup can provide topological defect acting as a seed, such a mechanism 
disallows possible eternal topological inflation. 

2.2 The modified WdW cosmology: results and observations 

In this paper, we use the general framework of large RG approach in matrix quantum me- 
chanics on a circle 12] that works with A^^ quantum mechanical degrees of freedom instead of 
TV eigenvalues and thus probes physics beyond the free fermionic approximation. The dilaton 
emerges in a more complicated way involving all the off-diagonal degrees of freedom and hence 
if is not simply truncated to ipo but is dynamical in all the modes. Using this approach we 
derive a generalized WdW constraint and the corresponding (all genus) wave function for FZZT 
branes with Dirichlet boundary condition on matter. Thus the quantum cosmology we study 
is that of one dimensional closed universe of arbitrary shape I = § cPae"^"^ with homogeneous 
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matter t = x and arbitrary scale fluctuation 51 ^ § £'ae^'^\pl5ip. 

In matrix quantum mechanics picture the universe of length / is created by the macroscopic 
loop operator jTr0^(t) that cuts a hole of lattice length / on the discretized string world sheet. 
They have the same continuum limit as the lattice loop operator iy(/,t) = -^Trexp(/ </)7v(t)/a), 
where </'7v(t) are N x N hermitian matrices and a is a lattice spacing, and can be studied in this 
form. The commonly used boundary condition on these one-dimensional universes is the one 
that constrains the boundary length to be / throughout which the matter field takes a single 
value (say t = x). It will be more useful to compute the wave function for the one point function 
of the loop operator with one puncture or the wave function for a boundary with one marked 
point 

^(/) = M , ^(/) = (^w{l, x)) = (AT-i Tr exp[/0(x)]) . (2.6) 

A large length fluctuation over an arbitrary shape of loops is the dominant effect at small scale 
and hence is crucial to explain the inflationary scenario of early universe. Below we summarize 
our results and critical observations. 

A nontrivial condensate 

The world sheet RG modifies the WdW constraint by adding derivatives with respect to 
scale of all degrees and introducing nontrivial matter dependence. Along the line of our dis- 
cussions above, this amounts to capturing physics of a nontrivial and complicated (perhaps 
nonintegrable) spatially deformed time dependent tachyon condensate. The inflaton field or 
(in this case) the matter field acquires a light (exponentially vanishing) mass and thus have 
no problem in executing a slow roll. The world sheet RG also adds a 0{1/N) correction that 
has inherently quantum mechanical origin coming from the anomaly in the Callan-Symanzik 
equation. 

The quantum creation of inflationary de Sitter vacuum 

The short distance / < H^^ ~ wave function, which we obtain, is a tunneling wave 
function describing quantum creation of a closed universe from nothing and its nucleation to an 
inflationary de Sitter vacuum 

m{l < H^^) ~ exp[-ci/2/2] . (2.7) 

At large scale / > H^^ the evolution of the universe is classical. Thus at large scale the wave 
function reproduces the semiclassical behavior of an expanding global de Sitter geometry, in 
spite of higher derivative and other quantum effects. The relevance of the higher derivatives at 
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short distance can be realized from the fact that, in this domain a quadratic WdW will give 
rise to a far past global de Sitter ( Jo(-\/— /i/), /i < 0) only, instead of a tunneling wave function. 




Suppression of baby universes 

The baby universes generated {i.e. emitted or absorbed by loop splitting and reconnection) 
in this cosmology are of vanishing sizes and their formation is suppressed by higher orders in 
Their amplitude comes from nonlinear part of the WdW equation, generated by multi- 
trace contributions in the matrix quantum mechanics loop dynamics, that can be represented 
schematically by a complicated '*' operation 



where H is the modified WdW operator as captured by the world sheet RG. The one point 
function ip is a Hartle-Hawking wave function in the sense of The amplitude corresponding 
to the operation is coming from some anomaly, breaking the scale invariance imposed by 
a linear WdW constraint Titp = 0, originating from nonlocal contribution from multiple world 
sheets. The vanishing size of the baby universes and suppression of their amplitude helps 
in continuation of inflation once it onsets. Note that, the birth of baby universes involves 
topology changing and thus inherently is a beyond minisuperspace process. Such amplitudes 
get contributions from the points in the interior of different topological sectors of the superspace 
We observe that the sub- leading contribution from the emission of a single baby universe 
self-tunes the de Sitter entropy to be small positive. 

Emergence of Hubble scale 

In 2D the gravitational coupling is dimensionless and there is no proper Planck scale. The 
Planck scale assumed to be spontaneously induced at ~ 1. However, in the large world 
sheet RG, the Planck scale emerges naturally as the Hubble scale for the problem. The potential 
in the generalized WdW equation for the one point function of loop operator with a marked 
point helps to identify the required inverse Hubble scale H~^, which is very small for c = 1 case. 
This would be favorable for quantum tunneling of universe from nothing through an inflationary 
scenario. The potential below this scale has a tunneling type barrier. The geometry at H~^, at 
which the universe pops out from the classically forbidden region, is a de Sitter geometry. 

The open string moduli dependence and the geometry of the true vacuum 

In our approach, the nonperturbative effects due to the A^^ quantum mechanical degrees 
of freedom is observed to be introducing explicit open string moduli dependence in the wave 




Hip + t(j * ijj = , 
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function^ and in the Hubble scale. Such a dependence is not present in the wave function when 
we explicitly work with quantum mechanical degrees of freedom. The open string moduli 
dependence does not directly affect the inflationary scenario, as the tunneling wave function 
solution exists for any sensible constant H . However, the explicit dependence of the Hubble 
constant on the open string moduli determines the geometry of the true vacuum one tunnels to. 
It will be extremely interesting to investigate whether the world sheet vortices are responsible 
for some kind of spontaneous symmetry breaking that chooses the ground state. 

Choice of initial condition 

A problem of describing creation of universe, in particular an inflationary universe, always 
has subtleties with choice of unique initial condition corresponding to 'nothing' ^Hl- This has 
long been a subject of big debate (see for example [H]). With introduction of more degrees of 
freedom as we go beyond minisuperspace and derivatives of all degrees, the situation is even 
worse. One would encounter an infinite dimensional initial condition (for example, specifying 
initial value of all the derivatives of the wave function) in all the degrees of freedom involved 
(/q, 51)- However we observe that a gaussian profile of initial scale fluctuation il){5T) ~ e"*^'^ pro- 
vides all the required initial conditions and thus uniquely determines a tunneling wave function 
at UV. In other words, the wave function for initial scale fluctuation can be viewed as a 

co-moving (gravitational) wave detector to capture the inflationary vacuum. 

The arbitrary scale fluctuations 

The arbitrary scale fluctuations, taking us beyond the minisuperspace, are playing a pivotal 
role in defining a proper cosmological picture here. Because of the scale fluctuation, the actual 
volume of the one dimensional universe can no longer be definitely small at weak coupling and 
large at strong string coupling. A large fluctuation of the form 5ip ~ ^^^v _ g+v^</' would drive 
the scenario to the opposite. Now the length scale of the universe / = § e'^°+'^ ^/2<p_g+^/2■p ]^g(,Qj^gg 
small as (/? — s> oo (strong coupling) and grows as ^9 ^ — oo (weak coupling). 

Amusingly, such a fluctuation in ip would imply a scale dependence of the power spec- 
trum of primordial density fluctuation for our one dimensional universe (naively, the density 
perturbation ^ ~ —50 ~ q+^v ^ ly Strikingly, the tunneling wave function fixed by the ini- 
tial gaussian scale fluctuation profile captures a vanishingly small power-like scale dependence 
(like a real- world nearly non flat CMB anisotropy). Let us define the probability of finding 
the universe between a scale I and / + 5/ as / > at a time 0+ after the big bang to be 
J = P{1 -^1 + 51) = ^'^li^i'l''^ ~ l"^^) ' ^ small universe (/ < 1) and for an estimate 
*We thank Micha Bcrkooz for a discussion on this issue. 
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/ ~ H ^'^^ to realize the bound I ^ H ^, the power spectrum is shghtly red-shifted 

n,~0.96. (2.9) 

P 

Arbitrary matter fluctuations 

An arbitrary matter fluctuation is crucial to capture the particle production and also the 
formation of large scale inhomogeneities in inflation (as seen by comoving particle detectors). 
In matrix quantum mechanics this is a more delicate issue as it would require to consider 
wave function for a Dl brane with Neumann boundary condition on both Liouville and matter 
coordinates. In some sense this implies dominance of the kinetic term in the matrix quantum 
mechanics action (hence presumably the dominance of the nonsinglet sector), over the potential 
term, that transforms as singlet under SU{N) rotation. The angular (off-diagonal) degrees of 
freedom in the kinetic term can not be integrated out except for large compactification radius 
or for matrix quantum mechanics on a line [HZlEHlEni- In other words, it would need a full open 
string description of matrix quantum mechanics. In [HZI the problem have been dealt to some 
extent by considering interacting spin chains. In an upcoming paper we will extend the 
case studied here to an inhomogeneous universe (the Neumann boundary condition on matter) 
to study arbitrary matter fluctuations. The case with Neumann boundary conditions in both 
matter and Liouville directions has not yet been solved in matrix quantum mechanics. 

The factor ordering ambiguity 

The arbitrary scale fluctuation allows a small universe at strong coupling. As / — 0, the 
potential in the WdW equation V{l,x) — > and hence the initial wave function could be 
incoming or outgoing plane waves (FRW geometries). However, large RG captures a factor 
ordering ambiguity p 1 of the quadratic differential operator {l^§il^§i) in the WdW equation. 
This bars the wave function from being only plane wave at small scale, thus encouraging an 
inflationary scenario. 

Instability and third quantization 

Intuitively one can imagine for eternal inflation, in which the space-time of universe nucleates 
in a fractal structure, a violent fluctuation in all modes of if would be necessary. On the other 
hand, such extra dynamical fields beyond minisuperspace usually give rise to unstable scalars and 
hence negative norm states leading to recoUapse of the universe jS] . However, in our approach, 
the situation is quite different. The theory of loops on world sheet can be seen as some kind 
of third quantization (see for example |^ 133] and references therein). Then the problem of 
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negative norm state is not relevant as the wave function of the universe is actually a quantum 
field operator (the loop operator itself) composed of creation and annihilation operators for the 
universe and the WdW equation follows from the third quantized action. The generalized WdW 
operator in large N RG emerges as a natural gauge condition for loops of arbitrary shape with 
arbitrary scale fluctuations (from the anomaly in the Callan-Symanzik equation for the loop 
operator), that takes care of negative norm states of the dynamical fields. 

The minisuperspace limit 

To discuss the validity of our derivations we reproduce the minisuperspace AdS2 quantum 
mechanics by large RG on the space of matrix eigenvalues. It is very interesting to note that 
the world sheet RG somehow captures a self-tuned deformation of the AdS2 quantum mechanics 
of the degrees of freedom (the matrix eigenvalues) to that of a ^5*2 quantum mechanics of all 
the de grees of freedom. Perhaps it is an indication that there is such a self-tuning mechanism 
of dynamical vacuum selection in real world that we will be able to capture by going beyond 
the minisuperspace approximation. 

D-instantonic picture for tunneling? 

It will be certainly interesting to understand if there is any instantonic picture for the 
tunneling. In free fermionic representation, the vacuum with de Sitter sign of the cosmological 
constant (/i < 0) has an instability due to an instantonic tunneling amplitude 0(e~^) from 
the other (AdS) side and is central to the nonperturbative formulation of the theory [HH]- The 
Liouville theory coupled to c = 1 matter does have D-instantons (ZZ branes) of c = 1 noncritical 
string theory localized in the strong coupling region {ip 0) [Hnj- It has been interpreted that 
the FZZT brane (extended in (/? ^ — oo) probes the ip direction up to a certain scale (until 
Qs ^ 1 analogous to / ~ in our cosmology), beyond which it dissolves. Below this scale 
the open string degrees of freedom are only these D-instantons stuck aX ip ^ oo which can 
be thought of as the remnant from the superposition of a FZZT brane stretching from weak 
to strong coupling and then back to weak coupling region [ZD]- Thus the wave function for 
FZZT branes at small scale / < actually computes a D-instanton tunneling amplitude in 
all genus. Note that, unlike the tunneling event in A^ degrees of freedom picture (eigenvalue 
representation), here the tunneling event captured by the dynamics of A^^ degrees of freedom 
ends in a stable de Sitter vacuum. 

An artifact of the formalism? 

One might wonder whether all these nice aspects captured by large A^ RG is in some sense 
an artifact of the formalism itself or they are really happening. To this end, we need to clarify 
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whether the small matrix field (small $Ar(t)) approximation, which is the only approxima- 
tion in the large RG, is in contradiction with the large eigenvalue at the strong coupling 
(eigenvalue D-instanton) or at the point of spilling ^, or more generally to a large quantum 
fiuctuation in the matrix field variable itself. However, in integrating out each momentum shell 
(^^ixAf , ^^'tvxi 5 = 0) from a (A^ + 1) x (A^+ 1) matrix $Ar+i(t), the remaining $Ar(t) part that 
enters in the interaction vertex $Ar is never approximated to be small in an absolute sense. 
Rather it is assumed lighter compared to the integrated out degrees of freedom which become 
heavier by a rescaling v v = y/Nv. This makes the interaction vertex lighter compared to 
external integrated out momenta, thus facilitating in the expansion of its exponential to get 
Feynman diagrams. 

World sheet vs space-time observer? 

Our study of the cosmology is performed with respect to the world sheet observer, which is 
meaningful as the observer should be a part of the universe. However, even though the wave 
function for the universe we are describing is computed on the world sheet, the loop operators do 
have a proper space time interpretation and the cosmology could be translated on that language. 

Euclidean vs Lorentzian picture ? 

It appears that the world sheet large RG with A^^ quantum mechanical degrees of freedom 
automatically captures the Lorentzian de Sitter vacuum of the supercritical theory instead of 
the Euclidean AdS vacuum of the subcritical theory, the later being captured by the large RG 
of A^ quantum mechanical degrees of freedom. Perhaps the world sheet RG has some nontrivial 
dynamics that leads to a growth of space-time dimensions^. 

In figure |21 below, we summarize the main results in a schematic form. The broken arrow 
represents a possible Woo deformations showing past and future spacehke infinities as was also 
included in figure HJ 

3 The large N RG in MQM 

We will now review in some detail the steps of the large A^ world sheet RG analysis in MQM. 
For a complete description see the original works in |T[|2j. The concept arises from the inter- 
pretation of the very existence of the double scaling limit (or the continuum limit) as some kind 
of Wilsonian RG fiow jHSj . 

^We thank Emil Martinec for pointing this out. 
^We thank Eva Silvcrstcin for a discussion on this. 
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Figure 2: Summary of results and observations. 

In the double scaling limit, as the matrix coupling constant g g^, the average number of 
triangles in triangulations at any genus G diverges as 

(nG)~(l-G)(7o-2)(l-(?/(?,)-^ (3.1) 
Simultaneously with N oo, the length of the triangles (the regularized spacing of the random 

1 2 

lattice) a N ^-70 scales to zero to keep the physical area a {uq) ~ 2-70(1 — g/gc) or 
equivalently the string coupling gs = N'^{g — gc)'^~^° fixed. The existence of double scaling 
limit indicates that a change in the length scale induces flow in the coupling constants of the 
theory in a way that one reaches the continuum limit with correct scaling laws and the critical 
exponents at the nontrivial IR fixed point determined by the flow equations. The large world 
sheet RG in MQM ^ |2] is basically the evolution of the two sets of parameters of the theory, 
the size of the matrices and the cosmological constant (mapped into the matrix coupling 
g) and all other matrix couplings, at the constant long distance physics with the rescaling of 
the regularization length in the triangulation of the world-sheet. In the Wilsonian sense this is 
done by changing A^ ^ A^ + 5N by integrating out some of the matrix elements, which is like 
integrating over the momentum shell K — dK < \p\ < A, and compensating it by enlarging the 
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space of the coupling constants g ^ g + 5g. Here the space of couphng constants will contain 
both the matrix coupling g and the mass parameter M^. 

3.1 The RG scheme 

Let us now schematically explain the emergence of the flow equations. By integrating out a 
column and a row of an (A^ + 1) x (A^ + 1) matrix, the c = 1 matrix partition function satisfies 
a discrete relation 

Zr,+i{g, M, R) = [\{g, M, R)f"z^{g', M', R') (3.2) 
with the following flow equations for the couplings 



9' = 9+j^PA9,M,R) + 0[l^) 



M'2 = M' + ^/3M^ig,M,R) + 0(^^y (3.3) 
and an auxiliary flow equation 

A(^, M, i?) = 1 + iV"^ r{g, M, R) + 0{N~'^) . (3.4) 

Here g',M' represent the renormalized couplings. Identifying as the world sheet scale of 
the RG transformation, the functions (3g,(3M are interpreted to be the beta functions for the 
corresponding couplings. The function A gives the change in the world sheet free energy (the 
string partition function) 

^(iV, g, M,R) = ^ In Z^ig, M, R) . (3.5) 

The auxiliary flow relation for (3x is useful in studying thermodynamical consequences. As 
A^ is taken to be large, the world sheet scale 1/A^ becomes infinitesimal. Then the discrete 
relation ()3.2|) for the Matrix partition function actually translates into a continuous Callan- 
Symanzik equation satisfied by the world sheet free energy with the function r controlling it's 
inhomogeneous part 

N— - (3{g, M,R)—- (3{g, M,R)-^ + 7] H^, g, M, R) = r{g, M, R) . (3.6) 

Here 7 acts as the anomalous dimension. We will see that the (3g,(3M computed by the large 
A^ RG are such that the homogenous part of the Callan-Symanzik equation indeed determines 
the correct scaling exponents for c = 1 matrix model around the nontrivial fixed point. The 
inhomogeneous part is related to some subtleties in the theory, like the logarithmic scaling 
violation of the c = 1 matrix model. 
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3.2 Integrating over the vectors (bosonic quarks) 

To carry out the explicit derivation of the flow equations, let us consider the usual hermitian 
matrix quantum mechanics for (A^ + 1) x (A^ + 1) matrices 07v+i(t) with a periodic boundary 
condition on the matrices. It has a cubic potential acting as a wall stabilizing the inverted 
oscillator potential. Now considering the following parametrization (the scalar being of relative 
order is neglected in the double scaling limit), 

, ... [ ^N{t) Va{t)\ 

and considering the power expansion (the higher order terms in v*v are suppressed by powers 
ofO(l/iV)), 

Tr0^'+i = TT(f)% + 2k v*(P'^-\^ + 0{v^) , 

Ti(t)''^Xl = Tr0^ + {2k + 1) v^cp'^-^ + 0(i;V''"') , (3-8) 

the matrix partition function reduces to a part containing N x N matrices and an integral over 
the vectors 



ZN+,[g,M,R] = [ V^'Mt) 

J</.jv{27rfl)=</.jv(0) 

exp[-(iV+ l)Tr^ ^ dt {^^^(t) + ^M20^(t) - |0^(t)}] 
/ V^v{t)V''v*{t) 

Jv,v*{2tvR)=v,v*{0) 

p2TrR 

exp[-(iV + l) / dtv*{t){-d^ + M^-g(l)N{t)}v{t)]. 
Jo 



(3.9) 



This is nothing but the c = 1 partition function suitable for Veneziano type QCD considered 
in jnni EZj to study open strings. Both color and fermion (quark) flavor A'^^ are taken to 
be large. Here the quarks are bosonic (the vectors) and A^ = A'^^. It is precisely these fields 
in the fundamental representation of the global SU{N) group, which generate boundary terms 
in the Feynman diagrams. Integrating adiabatically over these (bosonic) quark loops generates 
a logarithmic interaction term (+A^ Trlog[— + — g(j)fyr{t)]) and hence the determinant 
det|| - + - g(f)N\\~^ that gives rise to the boundaries in the world-sheet. Logarithm with 
minus sign would arise if the integration is performed on A^ flavors of fermions. Unlike the 
open string models mentioned above, here the couplings in front of the logarithm and in its 



20 



argument are nor introduced by hand, rather they are determined by the original closed string 
action. One can ignore the derivative term inside the logarithm if the mass and the couplings are 
large enough. This makes the dynamical loops uncorrelated in time and gives rise to Dirichlet 
boundaries on the world sheet. On the other hand, if the derivative term dominates then the 
world sheet boundaries are Neumann. In ^7\, c = 1 model with explicit expression for the 
fundamental fields has been considered and is the one which is closest to our model obtained 
after integrating out one row and one column of flavor degrees of freedom. If one considers 
infinite line, only ground state is relevant. For some choices of the coupling and the mass of the 
particle moving around the boundary of the holes, it has been possible to find the ground state 
jni], but the exact spectrum is not known. Our ultimate goal is to derive the exact spectrum 
of these boundaries and to study their cosmological implications. 



3.3 Evaluating det| | - df + - 




Let us now evaluate the determinant by standard Feynman expansion. In general it is a nonlocal 
object. For simplicity one could consider the constant 0Ar mode, which is equivalent to studying 
the effective potential to determine the phase structure. But here we will treat more general case 
of evaluating the determinant with flavors coupling to general (pj^j. The complicated induced 
interactions arising from the logarithmic term are ignored in small field approximation in order 
to get back the beta functions that depend only on the original couplings. In terms of the 
Fourier modes, the ^-integration can be expressed as 



1 f TTi^ 

I[g, M, R] = —=====— / ( n dvldv:) exp[- ^ v*^{'^+M^)v^+g v*^<Pm^iVi] . 

L V V /J n m m,l 

(3.10) 

The integral can be re-casted in terms of a gaussian part Iq being acted upon by a part generating 
interaction vertices 



I[g, M, <p^, R, N] = ^^^^ ^ exp [ - ^ 0^- J^?' <P)] ^o{R, M, N) , (3.11) 

mi,li 

with the gaussian part defined as 

/o [R, M, iV] = /" ( n dv*dv,) exp [ - 5^ v*^0::y^] . (3. 12) 

j m 

Here the operators O'^^ and O^Zi{g,4>) respectively define the inverse of propagator carrying 
a momentum mn/RSmn and interaction vertex carrying a net momentum {m — l)/R 
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/mn 



(3.13) 



The inverse of these operators define various propagators and vertices according to figure El 



m 
(a) 



m 



ml 



(b) 



Figure 3: The propagators and vertices: (a) [O] 



mnl ' 



(b) K^l.ig, 



The gaussian integration is basically absorbed in a net prefactor in ()3.11|) such that it 
reduces to the sum of one loop Feynman diagrams T,[g, M, (f)^, R, N] up to this prefactor. The 
sum Ti[g, M, (pN, R, N] is given by 






Figure 4: The diagrams contributing to S at one-loop order. 
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+ 0{g') (3.14) 



Now the one loop correction to the various coefficients in the effective action is evaluated term 
by term by inverse transforming the Fourier modes, summing up the set of infinite series and 
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then exponentiating and log-expanding T,[g, M, (pN, R, N] using small field approximation. The 
summation is done according to 



oo 



exp[i{m/R)t] _ ttR cosh^nM R - Mt) 



^ + M 2 M sinh ttMR 



0<t< 2tiR . 



(3.15) 



m=— oo 



Since after performing the inverse Fourier transform, the various terms has nonlocal integrals 
over several one dimensional dummy time variables, we breakup the variables into center of mass 
and relative coordinates. Then we expand the functions about the center of mass coordinates, 
assuming the relative coordinates to be small enough, and consider integration over the relative 
coordinates. The higher derivative terms are dropped due to smallness of the relative coordinate 
compared to the center of mass coordinate. For example, by the coordinate transformation 



(^1,^2) ^{T = (ti+t2)/2,r = (ti-t2)/2), the O^g"^) term of the form / dti j dt2 Tr[(j){ti)(j){t2)] 



comes from O(r^) term in the expansion, whereas the next higher derivative term appears at 
O(r^). All the mixed terms of powers of (p and the powers of its derivatives appear in odd orders 
of r and hence are automatically dropped due to integration of r over a symmetric range. 
Collecting all the contributions up to 0(000), the expression for I][g, M, (p^, R, N] becomes 



E[^, M, R,N] ~ 1 + Fgi{R, M) g dt Tr0^(t) + Fg^iR, M) g" \ dt Tr0^(t) 



Note that, in the evaluation of S we keep the contribution from the nonlocal terms of the action 
up to the kinetic term and the contribution from the higher order terms in the matrix field up 
to the cubic term. All other operators are redundant for our purpose and are negligible due to 
the small field approximation. The coefficients F{R)s are defined as follows 



can be expanded around T, considering r small (which is true for ^2 > 0). Then the 0^ term 
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Fg^{R,M) = -LcothvrMi?, 

FM, M) = \ i-nMR cosh 2tiMR + - sinh AtiMR) , 

1 /ttMR n'^M^R^ 

Fa2 (R, M) = ^ cosh An MR cosh 2n MR 

' M^ smh^n MR\ 16 6 

- — (1 + Sn^M^R^) sinh An Mr) , 

FaziR.M) = . 

^ 64M5sinh3 7rMi?(cosh27rMi? + cosh47rMi?) 

[AnMR{3 cosh nMR + 2 cosh 3n MR + 2 cosh 5nMR + cosh 7n MR) 

+ sinh vrMi? + sinh 3nMR + sinh SvrM/? + 2 sinh TvrMi? + sinh 9nMR] . 

(3.17) 

The tadpole term hnear in 0Ar(t) is now removed by shifting the background and setting it's net 
coefficient to zero. This fixes the value of the shift ( / ~ —gM'^Fgi/N + 0{^)) and accordingly 
modifies the coefficients of all the terms in the action. 

3.4 Restoring the original cut-off and the space-time geometry around 
a fixed point 

As usual in Wilsonian RG we now restore the original cut-off by rescaling the variables 



{t) ^ p(j)'^{t') , t'^t{l-hdL), R' ^ R{1 - h dL) . (3.18) 



Here dL = 1/N and /i is a function of the radius and the couplings whose functional form can 
be determined from the behavior of the flow equations near the flxed points. As one approaches 
a fixed point, h saturates to a constant characteristic to that particular fixed point, in the sense 
that it determines the scaling exponents. Thus the rescaling of the radius R in some sense 
tells us that as the system fiows to various fixed points, the target space geometry changes 
accordingly. This is something new in matrix model which directly enables us to determine the 
target space metric around a fixed point by solving the equation in it's neighborhood 

dR , „ , , 

— = -hR. 3.19 
dL 

For example, later in this section we will see that h = corresponds to a c = 1 fixed point. 
This trivial rescaling indeed gives a fiat metric {R^ = const). In [Z2], this principle is used to 
extract 2D Euclidean black hole metric from matrix quantum mechanics. 
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3.5 Flow equations and the nontrivial fixed point 

We now proceed towards determining the flow equations by setting the overall coefficient of the 
kinetic term to one. This gives p as 

p = 1 + ^(/i - 1 + g^Fg2)dL + 0{dL'') . (3.20) 

The resulting effective action is of the same form as the bare one with the renormalized strength 
of the couplings. The renormalized partition function is given by 



Z^^, = A' ^ Y '^^^^Nit') exp [ - iV Tr dt' + M^^^ - 

g' = g+[h-^^gdl+[Fg-s{R,M)-^Fg2{R,M)]gUL, 
M'^ = + [2hM^ + g^{l- M'^)Fg2{R,M) - g^M^Fg2iR,M)]dL. (3.21 

The resulting beta function equations are 



= %={lh-l)g+[FgsiR,M)-^Fg2iR,M)]g', 

Pm^ = ^=2hM' + g'[{l~M')Fg2{R,M)-M'Fg2{R,M)]. (3.22) 

As we already mentioned much of the structure of the beta functions depend on understanding 
the quantity h. The gaussian fixed point A* = (O, O) satisfies both the equations Pg = and 
Pa'P = trivially for any h. One can show that for /i = 0, or in other words for trivial rescaling 
(t' = t, R' = R) of the target space coordinates, the nontrivial fixed points A^i^g*'^ ^ 0, M* ^ 0) 
have c = 1 critical exponents. Such a pair of nontrivial c = 1 fixed points of the flow equations 
()3.22|1 exists for any R and are given by 



2Fg^{R,M*)-'iFg2{R,M*)' 
(1 - M*^)Fg2iR, M*) - M*^Fg2{R, M*) = . (3.23) 

The second equation solves for the values of M* for different R. The fixed point values of 
the couplings are scheme dependent. As such their explicit values do not affect the physics as 
long as they are finite. However, to study the Dirichlet boundaries on the world sheet, we are 
interested in large values of M*. As one cuts a hole on the triangulated world sheet, there is no 
longer any matter or the matrix valued fields sitting inside the cut portion. Instead there are 
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vectors v, v* who sew the cut portion by correlating the matter at the different points on the 
boundary. As their propagator goes hke 1/M*, a large M* would imply uncorrelated matter on 
the boundary and hence will serve as a Dirichlet boundary condition [HZj. For h = 0, solving 
the corresponding equation for M* fl3.23|) . we find that a large M* is feasible for irR < 1. This 
implies that the effects of nonsinglet sector cannot be suppressed in this situation. In fact M* 
can have any value at vanishing R and is of order one in the range < irR < 1. Around nR ~ 1 
it goes to a peak (M* ~ 2) and then rapidly goes to zero. 



3.6 The c = 1 critical exponent 

Let us now explain how a nontrivial fixed point of ()3.22j) with h = corresponds to a c = 1 
fixed point. The critical exponents for the scaling variables, the renormalized bulk cosmological 
constant A = 1 — g/g* and the renormalized mass (related to the renormalized boundary 
cosmological constant) Ai = 1 — M/M* can be determined from the eigenvalues of the matrix 

The homogeneous part of the Callan-Symanzik equation, satisfied by the singular part of the 
free energy, can now be rewritten as 

The eigenvalues fijS of the matrix fi^ / are nothing but the scaling dimensions of the relevant 
operators. The singular behavior with respect to the renormalized bulk cosmological constant 
goes as, 

J^s ~ A^/^Vi[iV^^ A ] f2[N^^ M ] . (3.26) 
Comparing the above expression of Ts with the matrix model resuh ~ A^^-'ro) /[Ar2/7iA], 
or using the standard definition of the susceptibility V ~ ^^^1 M=o ~ A^'^^ , the string 
susceptibility exponent 70 is given by 

70 ~ (2 - 2/fii) . (3.27) 

Note that in our analysis 71 ~ 2/f2i is consistent with the matrix model relation 70 + 71 = 2. 
This relation is independent of the explicit values of 70 and 71 and is easily obtainable from the 
consideration of the torus. The string susceptibility exponent at genus G is defined by 

7G = 70 + G 71 . (3.28) 

For the nontrivial fixed point, Vtu = 1 + 5/i/2. Note that in this RG, typically g* is small. 
The pair of nontrivial fixed points is always situated close to the Gaussian fixed point. This 
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implies the off-diagonal elements Qu, ^21 , being proportional to the powers of g*, are small and 
Q22 ~ h. The scaling dimension matrix is effectively a diagonal one with Qi = l + 5/i/2, Q2 = h. 
Thus h = gives the c = 1 critical exponent, 70 = 0. This is true for any R. The pair of c = 1 
fixed points we thus get are repulsive with respect to the flow of the parameter g while the 
gaussian fixed point is attractive. One can show that this pair of fixed points reproduces all the 
known physics of the c = 1 theory and respects T-duality. 

4 The AdS minisuperspace from large N RG 

We would like to know that in the context of our large RG, which is formally dealing with 
the full superspace, is there really any limit in terms of the scale / that would confine us to the 
minisuperspace and reproduce the familiar minisuperspace results computed by the continuum 
Liouville theory or by the free fermionic matrix model? We would also like to know that if such 
a limit exists, then how the full quantum wave function is differing from the minisuperspace 
wave function away from this limit. This will help us in understanding the limitations of such 
a truncation of superspace in computing the actual quantum wave function of the universe. 
We will see that the minisuperspace physics in c = 1 matrix model corresponds to quantum 
mechanics of degrees of freedom, the eigenvalues. Here, we will show that the large RG 
acting on matrix quantum mechanics observables on the space of eigenvalues correctly captures 
the familiar minisuperspace WdW equation. This tests the correctness of the approach in 
deriving the WdW constraint. 

To be more specific, let us consider the one point function of loop operator in the space 
of + 1 eigenvalues from the diagonal matrices Ajy+i^x) = Q~^(f)N+i{x)fl, Q G SU{N) with 
periodic boundary condition A(0) = \{27iR), given by 

(j'j.qI-^n+i(x)^ = (^e''^^^^)^'^'^ . (4.1) 

We then integrate out one eigenvalue Ao(x) adjusting the matrix couplings according to double 
scaling limit, thus keeping g^ fixed. This generates a difference equation 

(j^e^^M+iix)^^^ _ (Tre''^^^^^^ = ((e'^«(^\^ - l)(Tre'^'^(^)) , (4.2) 
where, 77- are the renormalized matrix couplings given by 

g' = gil-l/2N), (3, = -g/2, 

M'2 = , /?M2 = . (4.3) 

Note that, the beta functions above are too trivial to capture the nontrivial c = 1 fixed point 
and thus only have the trivial gaussian fixed point. This happens because in the space of the 
eigenvalues one practically works with the gaussian matrix model. 
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Now the left hand side of the difference equation can be interpreted as the Callan-Symanzik 
operator 



( W^--(^)),, - (Tre'^-(^)),, = - Sg-^) {Tre^^^^'^ , (4.4) 

that simphfies to 



= n(- f dt[X{t)l-{-d^ + M^)X{t) - gX%t)y^^^A (e'^^'M-^ . (4.5) 

\ J 2 I ri'. 

Now, just hke in free fermionic theory, rescahng Xn ^ cubic interaction term becomes 



0{1/N ) and we focus on the top of the quadratic potentiaL Simultaneously we rescale I ^ y Nl 
keeping the arguments of the exponentials unchanged. This makes 

" (Tre^^"^'^).^ = -\M'd^{Tre^^-^%. . (4.6) 

Plugging this in (j4.2|) we have the minisuperspace WdW equation 

[~MH'^f + A^^l']m=0. (4.7) 

Here the wave function ip{l) is given by the one point function (Tre''^^^^-')^'. Comparing ()4.7|) 
with the standard minisuperspace WDW equation (|2.4p . the world sheet cosmological constant 
is given by 

4/i = 8/x/M2 = A(i_(e'^o)^,). (4.8) 

Now to determine the sign of the cosmological constant, recall |13 |^ that in the free 
fermionic theory the single particle wave functions in the inverted harmonic oscillator potential 
are supported at — oo < A < (note that, here we have opposite sign of the cubic interaction 
term as that of jlH])- Also the cosmological constant being of 0(1/M^) is small. Thus /i ^ 
(/i ^ as Ao ^ and g ^ 0) and comparing ()4.7|) with ()2.4|) . we see that the large N RG 
of the matrix model observables defined on the space of eigenvalues indeed gives rise to AdS 
quantum mechanics, as would be expected from free fermionic theory. This tests the validity of 
our approach for the WdW cosmology we are addressing in this paper. 

5 The modified Wheeler-de Witt constraint 

So far we have discussed the setup for large RG analysis giving rise to flow equations that 
have nontrivial c = 1 fixed points. We also have discussed the minisuperspace limit of the large 
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RG flow that gives rise to the famihar AdS QM computed by free fermion representation 
of the c = 1 matrix model. Let us now move towards the computation of the wave function 
by deriving the modifled Wheeler-de Witt (WdW) equation that dynamically determines the 
required wave function. In fact, the Callan-Symanzik equation for the one point function of the 
loop operator with Dirichlet boundary condition itself gives rise to this constraint. 



5.1 Integrating over the vectors 

Let us consider the expectation value of the loop operator at a flxed x with respect to the 
hermitian matrix quantum mechanics path integral with a cubic potential and couplings rji = 

^(0 = ^(Tre'^(^)>^^. (5.1) 

As before, the matrices obey periodic boundary condition. Considering the parametrization 
(13. 7p and expanding the operator, we have 

DMt) e-~^^^^h,^t) Yl ^[Tr0^(^) + kv*^{£)^%-'ix)v,,{x)] . (5.2) 

fc=0 

Here Si\f{(f)) is the v,v* independent part of the reduced action given by 

I .12 T\/r2 

A 

The integral over the vectors I^[t) is the usual one as in ()3.9p 



S^{ct>) = N Tr(l + + ^<Pl - l<Pl ■ (5.3) 



/l'2nR 
dv*j^{t)dvN{t) exp[- J dt v*j^{t){-d^ + M^- gMt))vN{t)] ■ (5.4) 

As we have discussed already, this is basically similar to the integral over the bosonic quark 
loops in Veneziano type QCD with a large color and flavor (A^ = Nf). The change in the 
expectation value of the loop operator due to the renormalization of the couplings rji can now 
be expressed as 

(Tre'*-+^(^)),, - (Tre'<^-(^)),. = /(^ /,(t)t;^(x)e''^-("V^'(£)^Jv(x)) . (5.5) 

\ 2^ / ri'- 

Here S represents the diagrammatic expansion of that eventually renormalizes the bare 
action (recall (j3.16|) ). Thus the expectation value of any operator O with respect to the renor- 
malized couplings 7]'^ can formally be given by 



(O)^, = j D<i)N O e-^'^l*! S . (5.6) 
We will now closely examine the meaning of the both sides of the equation ()5.5|) . 
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5.2 The Callan-Symanzik operator and the modified WdW con- 
straint 



The change in the expectation value in ()5.5p is essentially the Callan Symanzik operator in 
discrete form, acting on the one point function of the loop operator. Consider the discrete 
relation 

(Tre'<^"+^(^')),, - = - H^] (Tre'<^"^'^),, • (5.7) 

Now recalling ()3.3|) defining the beta functions = N6r]i in the large RG, and using the 
identity J^(Tr e^'f'^^^^)^^ = {Nd/dN + l)^(Tr e''^'^^*))^^, the left hand side of dES)) can clearly 
be written as a Callan-Symanzik operator acting on tp{l,x) 



d ^ d 

However, the wave function is not completely annihilated by the Callan-Symanzik operator. 
The nonvanishing right hand side is essentially computing the anomaly in a nontrivial way. The 
origin of this anomaly lies in the fact that the loops of same length can continuously be mapped 
into each other by 2D diffeomorphism which acts like a gauge symmetry. The WdW operator 
acts as a gauge fixing condition. 

We will compute the anomaly below and will see that it is proportional to (j)]^^{x) which is 
singular in the small field approximation. The singularity can be tackled by taking derivatives 
with respect to / on both sides of the Callan-Symanzik equation. The right hand side then 
simplifies to a differential operator acting on the one point function. In this form, the operator, 
which annihilates the one point function of the loop operator, gives nothing but the modified 
WdW constraint. Physically this constraint governs the dynamics of loops of arbitrary shape 
with arbitrary length fluctuations as compared to the usual (minisuperspace) WdW constraint 
that only deals with the dynamics of a spherical loop. We will discuss in detail in what sense this 
constraint is 'modified', i.e. 'beyond minisuperspace', compared to the usual minisuperspace 
WdW constraint and it's implications. In this paper, we will show that the very existence of 
such a general 'beyond minisuperspace' constraint in low dimensional noncritical string theory 
is the key to probe question like inflation and topology change. 



5.3 Computing the anomaly 

Let us now go back to the computation of the anomaly term. As before, let us rescale the 
vectors as Vi — > Vi/V^nR and Fourier decompose all the vector and matrix degrees of freedom. 
We now perform the following small field expansion of the non-gaussian part of ()5.4|1 containing 
the interaction vertex = g(t>m^i 
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m,n m,n p,q 



.(j)p-qVq + . . . ) exp[- ^ v*^{m^lR^ + M^)v, 



'm 



m 



We then contract all possible v* pairs in the expansion that (upon acting on the remaining 
gaussian part of the integration) extracts inverse propagator (Cmn)~^ — + ^^)~^<^m,n.- The 
gaussian integration /q gives an overall pre-factor that neutralizes a part of the denominator. 
The contractions nontrivially mixes up the momentum modes at interaction vertices generating 
complicated interaction terms in the diagrammatic expansion. However, the Fourier modes of 
the singular term e^'l'^'^^^(f)'J^{x) do not get mixed with that of other terms. Hence the term is 
effectively decoupled from the rest of the expansion and does not receive any correction. The 
rest of the term with quantum corrections can finally be cast as a series summation of differential 
operators which we will explain below. This decoupling property of the singular term eventually 
gives back the one point function of the loop operator and is thus crucial in closing the WdW 
constraint. Here for simplicity we do not show the Fourier decomposition of the term. The 
resulting diagrammatic expansion due to the integration over vectors with insertion of the loop 
operator at the Dirichlet boundary is 





,i(m—n)x/ R 



m—n 




n—q 



_ ^{m—q)x/R 



9<P: 



m—n 



Figure 5: Diagrams representing 0{g) and 0{g^) terms in ()5.10p . 
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+ ... 



(5.10) 



where S denotes the non-Gaussian part of the diagrammatic expansion of the vector integral 

with /q representing the Gaussian part. It is now convenient to inverse Fourier transform the 
expansion term by term and then evaluate the Feynman diagrams by summing up the series 
with exponentials. Thus in terms of the target space diagrams, the small field expansion for the 
anomaly becomes 



^ ^iTrp<^^(^V^'(£){^coth7rMi? 



27ri?\E 



M^J 2nR^V^ Rm^J ^^'^V^Rm^J 

m n 
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M2/22 ■ ^ 



+ 



6M^y ImF^^I^+M^J ' '^^^i) 2.^2/^2 + ^2'^^^^ 
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Now to evaluate the Feynman diagrams, we use the open string moduh x to replace the inverse 
propagator by a differential operator of x 

This simplifies the summation over exponentials in the above expression to delta functions 
biti — x) that forces any matrix quantum mechanics time tj (matter) to be at x, which is 
basically reflecting the fact that we are talking about Dirichlet boundaries that are introducing 
explicit dependence on the open string moduli due to nonperturbative effects coming from all 
the A^^ quantum mechanical degrees of freedom. Note that replacing the inverse propagator by 
f)5.13|) is not only consistent formally but also can be seen to be in agreement as an expansion 
in the large M limit for Dirichlet boundaries. In the extreme M — > 00 limit, the explicit open 
string moduli dependence is washed out. After performing the delta function integrals over the 
matrix quantum mechanics time, the diagrammatic expansion can now be written as 



(Tre''^^+^(*))^^ - (Tre''^'^^^'))^^ 



/ / 1 



— [^ — cothvrMi? Tr(e''^^(*V^^(5)) 



+^Tr(e'<^-(^')0]^i(x) O.mOx) + ^cothvrMi? ^ Tr(e'^-(^)0]^i(x) 0,<f{x)0,) 
Tr(e'<^-(^V^'(ii)C?x0(x)C»£(/.(f )C»£(/.(x)0,) + ) , , (5.14) 



The above diagrammatic expansion can now be arranged according to even and odd powers 
of matrix couplings which sum up to even and odd series of -^Oj^(^[x\ thus computing the 
anomaly of the form 

^rj^g;<^^+i(x)^ - (Tre'-^-^")),. = ^{t^^ Tr[(e'*-("V}v'(£))('^i + '^2)]) , (5.15) 
where the diagrams with even powers in coupling are given by 



^1 = ^cothvrM/2[l + i^O,,V(x)^0, 

= ^cothvrMi? cosh(^O,0(x)), (5.16) 



and the odd powers in coupling are given by 
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(5^17) 



Note that, adding up the series neatly to hyperbohc functions is consistent with the small field 
approximation as well as the large M limit for the Dirichlet boundaries. 



^1 = 




+ 




+ . . . 




Figure 6: Diagrams contributing in 5"! and 82- 



5.4 The modified WdW constraint 

Let us now differentiate both sides of ()5.15|) with respect to the length / of the loop to remove 
the singular term from the right hand side. Also we do small field expansion of Let 
us now recast powers of (f)N{x) to be multi derivatives of / on the loop operator. This closes 
the WdW constraint by factorizing the right hand side of ()5.15|) into a differential operator of / 
and X operating on the one point function of the loop operator and a modified WdW equation 
emerges. Thus 



d_ 

m 



(-Tt e''^^+^^*^') — (-Tt e^'f''^^^^\ 

\] /ri, \] lr)[ 



2M 
1 



coth 71 MR cosh 



VM2 ^9/ 



+ _ sinh ( 



(Tre'*(^))^. + . . . , 
(5.18) 



where the ellipsis corresponds to the multi-trace or worm hole deformations that leads to topol- 
ogy changing amplitudes. These are redundant operators and are in relative 0{1/N) compared 
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to the leading contribution and contribute to topology changing transitions. Using ()5.8p the 
linear (dealing with loop of a single topology) generalized WdW equation can be written as 



1 

P 



+ 



ld_ 
Idl 




coth ttMR cosh 




M2 



Ox sinh 



(5.19) 



Clearly the WdW constraint is modified or generalized, as compared to the usual one, in the 
sense of being beyond minisuperspace (see the discussion in section 2). It captures effects of 
derivatives of all orders, which is more than the effect caught by a perturbative superspace 
(considering fluctuations on the Fermi surface in the free fermion language). Also the higher 
derivative effect captured here is for all genus wave function as compared to the genus zero 
result obtained in the perturbative superspace. Moreover, the scale factor / here contains all 
the modes of the metric fluctuation as compared to the presence of zero mode /q only in the 
dynamical Fermi sea picture. Note the difference in the / dependence of the coefficients of 
V{d/dlo) and V'{d/dlo) in ()2.5|) and that of the cosh and sinh terms in (j5.19|) . This point needs 
to be understood. Another aspect of ()5.19|) being beyond minisuperspace is that, it has an in 
built way of determining the factor ordering ambiguity. However, as we mentioned above, the 
linear modified WdW constraint that governs the dynamics of loops of a particular topology, 
comes from the leading order contribution of The sub-leading orders, beginning from 

a relative 0(1/A^) compared to the single trace terms, give rise to the nonlinear part of the 
WdW constraint corresponds to dynamics of topological transitions between loops of different 
topologies. This is strictly a beyond minisuperspace process that can no way be captured within 
minisuperspace. We will return to this in the last section. 

From now on let us write ()5.19|) in terms of the wave function (j5.1|) . In order to reduce the 
dimension of the space of independent variables, let us consider the following ansatz for the 
Callan-Symanzik equation of the macroscopic loop variables 



{N—-P, — + l)i;il;N,r|,,R) = {l^J^,{'ni,N,R)+J^2iv^,N,R)) ^P{1; N,rii, R) . (5.20) 



Here, [J-'i] = [/~^] and JF2 is dimensionless. The ansatz is chosen keeping in mind that the wave 
function behaves somewhat like the exponential of the Euclidean action, the area P. From 
our discussion in the previous subsections we already know that the macroscopic loop has an 
anomaly in the Callan-Symanzik equation which demands that they are annihilated by the WdW 
operator that acts as a gauge fixing condition for them. Here, by considering this ansatz, we are 
actually assuming that we are looking at those macroscopic loops which not only have WdW 
constraint as a gauge fixing condition, also they themselves are invariant under scale changes. 
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So that, for them the anomaly basically simplifies into some kind of anomalous dimension (a 
number). This reduces the space of independent variables to the dynamical variables only, the 
length of the loop /. As we mentioned before, due to nonperturbative effect coming from all 
the iV^ quantum mechanical degrees of freedom, the wave function has explicit dependence on 
the open string modulus x which acts as another dynamical variable. Now our job is to plug in 
fl5.20|) into ()5.19p and solve for the wave function that as well determines J^i and in a self 
consistent way. The functions J^i and ^2 keep the information about the c = 1 fixed point of 
the matrix quantum mechanics around which the solution is being studied. 

Let us consider another simplification. Since the time direction is compact, ct = can be 
identified with x = 2ttR. Therefore one can assume the x dependence of the wave function ■?/;(/, x) 
to be periodic. The simplest choice, that is nonzero at x = 0, is plane waves il){x) = cos(^x) 
with quantized momenta. Thus, 

= A„ ^(x) , K = f^,^^, . (5.21) 

Thus considering (j5.2Uj) . (j5.19j) and (j5.21|) . the beyond minisuperspace WdW equation for the 
one point function of loop operator becomes 

21 ylcosh (g A„^) + 21 EA^sinh (Q A„^) - (/^^i + J^s)^ - 2Z^i]^/'(/) = , (5.22) 
where, 

A^^.otb.MR, B = |j. Q = -|j. (5.23) 

It is important to note here that all the parameters A,B,Q have M — s> —M symmetry (see the 
discussion under c = 1 fixed point in the previous section). 

To study the solution of ()5.22j) at different scales, especially at small and large scale limits, 
it will be useful to change the dependent variable to one point function with one puncture 

ib(l x) 1 e''^^(^) 
n,£) = ^ = j^{T^r'-^),,. (5.24) 

This basically defines a closed Dirichlet boundary of length / with one marked point on it. Using 
the following identity 

e^«- ^ = 7 m + (5.25) 

in ()5.22|) . the WdW equation for the one point function with one puncture can be rewritten as 

21 A cosh (g A„|) +21 BK sinh {Q A„|) - {f:F, + r,)j^-2ir, - ^^'^^^^'g^A^^"^ ] *(^) = ^ " 

(5.26) 
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This actually induces a polynomial potential term in the WdW equation that brings in the 
notion of a special scale l{n) = QXn separating the UV and IR physics. From 2D cosmology 
point of view this is naturally a inverse Hubble scale for the problem. In the next section 
we will discuss the significance of this scale in detail and study physics of the solution to ()5.26|) 
with respect to it. 

6 The cosmological implications of the wave function 

In the previous section we have seen that the modified WdW equation ()5.26|) for one point 
function with one puncture has a polynomial potential that naturally brings in the notion of 
a special scale l{n) = QXn for the problem. In this section we will study the cosmological 
implication of this scale and the study the small and large scale behavior of the solution to 
fl5.26|) with respect to this scale. 

6.1 The inverse Hubble scale 

Let us consider a rough qualitative analysis of the shape of the polynomial potential in ()5.26|) 
with respect to the scale /. Here one can expand the hyperbolic functions of differential operator 
and write down the potential with respect to the quadratic differential term (— /^^) in the 
Liouville metric dip = dl/l. Such an expansion is possible considering Q = g/M"^ to be small, 
which is definitely the case around a c = 1 fixed point that has small g and for Dirichlet 
boundary that has large M. The potential is given by 



Here n is related to the momentum modes of matter. Since [J-'i] = [A] = [/"^], the ratio Ti/A 
is purely a number. Now the structure of the potential depends on the size of the quantity 
BXn/A. For Dirichlet boundaries 

BXn/A = 2Mtanh(7rMi?)/(M2 + _) ~ 2/M « 1 . (6.2) 

R-^ 

This shows that for / << /(n), the potential V{l,n) looks like a barrier in a tunneling problem. 

This is more clear from the form of the potential ()6.1|) along with the constraint ()6.2|) in x 
space 

/, X 9^ cosh JM^ - gyP(7rR-x) , , n / n 

V{l,x) = Viocai + ^^==== 2 ,, ^^(2 - ^TTtanhTTMRl/g) , (6.3) 

^y]\/P - g^/P smh a/M2 - g^/P ttR 
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where 



2/2 



Viocai = — [2 - {C-^ + 2Ce-^)][M^6{x) + 6"{x)] + (2 - 4M tanhn MRl/g) 6 (x) , (6.4) 



being proportional to 6{x) can be neglected by considering x ^ 0. Thus the magnitude of 
the matrix quantum mechanics parameters around the c = 1 fixed point as compared to the 
smallness or largeness of / and the form ofV = V — Viocai controls the behavior. Let us consider 
that, as / — > 0, (7// <^ M. Then, assuming x ~ vri? ~ 1, simplifies to 



On the other hand V{1 00) -00 and V ^ at a scale / 0{g/M) > I ~ 0{g/M^). This 
provides the tunneling barrier. Recall that the special scale / occurs before the potential goes 
to zero. 

This indicates that for small scale / < /, the ID closed universe with homogenous matter 
that we are considering, presumably admits a tunneling type of wave function that creates 
the universe from nothing by Euclidean tunneling. In cosmological models of 4D gravity such 
tunneling wave function have been shown to predict infiation Whether the tunneling actually 
predicts an infiationary scenario depends on the scale to which universe tunnels to. If the scale 
expands exponentially in a short time, then the universe undergoes infiation. At large scale 
/ > / the potential approaches —00 and hence the wave function of the universe presumably 
admits outgoing plane wave that gives an expanding de Sitter universe in the far future. In 
the following subsections we will show that the phase / < / will be the phase at which universe 
undergoes infiation and the phase / > / will be the far future phase of an expanding universe 
which is already hugely infiated. The region, where / = /, actually gives the hot big bang. 

Thus the scale l{n) acts like an inverse Hubble scale that separates the small and large 
scale behavior of the wave function. For c = 1 matrix model, it is at very short distance as would 
be appropriate for an infiationary scenario. Note that the nonperturbative effects due to the A^^ 
quantum mechanical degrees of freedom is introducing explicit open string moduli dependence 
of the Hubble scale. Though the moduli dependence does not directly affect the inflationary 
scenario, as one always has a tunneling wave function solution for any sensible constant H, 
its explicit dependence determines the geometry of the vacuum one tunnels to. The matter 
dependence of this scale can be written down by inverse Fourier transforming l{n) as 
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(6.5) 





(6.6) 
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This is the scale factor corresponding to a de Sitter geometry of the world sheet 

„2 d2 

ds"^ = dX^ - cosh2(7rMi? - MX) d(p'^ . (6.7) 

Recall that the world sheet matter X{a) is mapped to the matrix quantum mechanics time 
through the delta function in the Dirichlet boundary condition ()2.3p . The issue of the Lorentzian 
signature is related to the fact that the large RG with A^^ matrix quantum mechanics degrees 
of freedom captures the de Sitter sign of the cosmological constant (/i < 0) in the WdW equation 
(like the supercritical case) instead of the usual AdS sign (/i > 0) of the c = 1 quantum 
mechanics. This effectively implies an imaginary slope of linear dilaton (7^ < , 7 — > 27 , Q = 
7/2 —s> iQ) and hence indirectly implies growth of number of space-time dimensions. On the 
other hand the exponential potential wall e'^'^ remains the same in the WdW equation. Hence 
with 7 — i> Z7, the Liouville field also becomes imaginary ip — > —iip. This takes care of the reality 
of the Liouville action and gives rise to a Lorentzian metric. 



6.2 The small scale behavior: inflation 

Let us now solve for the small scale behavior of the wave function. For I < I = QXn, The 
polynomial term can be expanded as 

9/2 _ °° 72m 

P(l)^ (Al-BXJ)-£p;;- (6.8) 

m=0 

As indicated by the potential the wave function in the region / < / can have tunneling 

through the barrier. Thus an appropriate trial solution satisfying (j5.26p would be 

^(0 ~a(A^,r/,)e-'^(''^)'' (6.9) 

Later we will argue that such a choice of trial function is a unique choice fixed by the initial 
profile of the scale fluctuation ^'(O + SI). Now using the identities 

cosh([^)e-^'' = cosh(2A//)e-^(''+'') , 

sinh([^)e-^'' = sinh(2A//)e-^(^'+^'') , (6.10) 

and expanding the hyperbolic functions in power series of / (for I < 1 i.e. g < and 
for a small I), we set the coefficients of all the powers of / to be zero. Thus we see that 
\Ef(/) ~ a{N,r]i)e~'^''^^^^'^ is a solution of (j5.26|) with the spectrum A determined by 
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A A 

A 9 79 A/'a A^i 1 

/2 g-A; ^ — i ^ ^ ^ Q 
2/2 

(^^e-A^^ + (_l)"i = 0, n^0,2. (6.11) 
n! /" 

Here we have three equations and three unknowns A,jFi and 0^2- The equations ()6.1ip indi- 
cate that any solution for the spectrum A should go somewhat like LambertVT function. For 
Lambertiy functions the function itself is in the same footing as it's exponential. Hence such 
a solution for the spectrum A will be a good candidate to actually perform the cancelation of 
the hyperbolic function by a polynomial ()6.8|) needed to satisfy equation ()5.26|) . Thus solution 
to the spectrum A can be written as A = C//^, where the number C is given by the three sets 
of relations in combined form 

_c n!(-l)"+i ^ 1 , ^ J-i CJ^2 



Now solving these relations we have 



C = -n LambertW ^ — ), n ^ 0, 2 , 

V 2n 

^2 = AP (-^-^-2e-^)- (6-13) 

Thus, A is given by 

A = -^ Lambertw( ^~^||'^" ) [m^ + n^0,2. (6.14) 

In order to give an idea, the table below illustrates different values of C corresponding to 
sample arguments, using (|6.13p . 
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n 


C 


1 


-0.79+0.77i 


3 


1.5 


4 


0.62+1. 13i 


5 


1.91 


6 


1.4 +1.27i 


7 


2.37 


8 


1.49+0.8i 


9 


2.84 


10 


2.68+1. 26i 


11 


3.69 



However, from ()5.20|1 A, jFi,^2 and the amplitude of the wave function a satisfy additional 
constraints 



^ OA _ 
Pvi — -Fx, 

OT]i 

In writing down this constraints we considered the fact that all the four parameters A, JF^, JF2, a 
are independent of / and A does not have any explicit dependence as indicated by the 
solution (|6.14p . Obviously the second one of the two relations in (j6.15j) determines the unknown 
parameter a. Now we need to check whether the other one is really a redundant condition or 
a serious constraint on the already determined parameters. To check this, let us recall that 
around c = 1 fixed point we have /i = 0, i.e. j3g ~ —{g — g*)/2 and Pm'^ ~ 0. Thus, a.s g ^ g* 
around a c = 1 fixed point, solving (|6.15|) for the ansatz A = TiA[g) gives 

Here K is nothing but an undetermined integration constant. Choosing this constant to be 
K = C{M'^ + makes the solution for A in ()6.16|1 equivalent to that already given by 

()6.14|1 . Thus the first of the two relations in ()6.15|) is a redundant constraint for A. Now let us 
solve the second relation in ()6.15|) by separation of variables, a = jj^ a,{g, M, R), which (up to a 
constant factor) gives 



9l 



Mg2 



COth(TTMR) 



1 



-2e~ 



^)coth(7rAf_R.) 



1 



-2e- 



(6.17) 
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Here C is given by the spectrum ()6.13|) . We see that (in accordance with [321) the multiphcative 
prefactor a absorbs all the non-universal parameter dependence of the wave function and the 
universal part goes like < /) ~ e . Note that, the R dependence in the / comes only 
through the Fourier modes of matter indicating matter dependence. Thus the full solution for 
the wave function at short distance {i.e. below Hubble scale) can be summarized to be 

^{l<l) = j^{g-g*)^ " ^ 

xe^^ ^ J e-^'/', (6.18) 

where C is given by (j6.13|) . Note that the nucleation probability given by this wave function 
is a mixture of the usual exponential suppression and the power law /i^ (here /i ~ ((7 — g*)), 
exhibited by the inflating lower energy phase in the critical droplet fluctuation in 2D Liouville 
gravity, that tends to mimic the behavior of dynamical lattice Ising model j25 • 

The first few values of C computed in table above shows that for > 3 the wave function 
indeed is a decaying wave function. In this regime, for odd n, C > and the wave function 
is a purely decaying wave function. For even n, C is complex. Hence the wave function has 
a outgoing plane wave part along with a decaying part. Note that, only for n = 1, the wave 
function has a growing behavior (the usual Hartle-Hawking behavior) instead of a decaying one, 
along with an outgoing plane wave. However, n here are in one to one correspondence with 
the orders of the derivatives of /. For computations away from minisuperspace, indeed n > 3, 
and hence there are enough spectra giving a decaying solution. The growing solution can be 
suppressed by the choice of the boundary condition \I^(0) = \E'o 7^ 0. Another point to be noted 
here is that, since for any value of n the wave function does not have any incoming plane wave 
part, one can say that the closed universes described by our wave function never re-collapses. 

Such an exponentially decaying wave function can be interpreted to describe quantum cre- 
ation of a closed inflationary universe filled with homogenous field x tunneling from a state of 
vanishing scale or nothing. To understand what kind of vacuum it tunnels to, let us consider 
the functional dependence of / on x from ()6.6p . At the end of the tunneling, eventually at a 
scale l{x), the universe in a state of large total energy (for example note that, in the stan- 
dard chaotic inflation scenario, the initial energy density is of the order of M^) E M (here 
[M] = \p] = [x~^ = [E]) will undergo a huge exponential expansion over a very short period 
of MQM time M'^ or due to a very small quantum fluctuation 6t ^ M ^. This is indeed 
happening for Dirichlet boundaries {6x ~ 0{1/M) 0). At the Hubble scale, to which the 
universe tunnels to, a very small quantum fluctuation sets an exponential growth in scale for a 
universe of huge energy. Thus for I < I, \6.1^) is the quantum wave function representing an 
one dimensional homogenous universe tunneling from 'nothing' to a state of chaotic inflation. 
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Note that the Hartle-Hawking wave function computed away from the minisuperspace ap- 
proximation, predicts 'tunnehng to an inflationary vacuum from nothing', in a structurally sim- 
ilar way (apart from the fact that here / = /(0(r)) and not just / = /o(^)) to the Linde-Vilenkin 
tunneling wave function in minisuperspace cosmological model 

^(/o < H^^) ~ e-5^''o . (6.19) 

This is unlike the usual growing behavior of the Hartle-Hawking wave function that makes 
predicting inflation to be impossible. One can argue that the boundary condition \t^(0) = 7^ 
chosen here may introduce arbitrariness in the initial condition of the universe. However, the 
Hartle-Hawking wave function '^{l) is the wave function for the one point function with one 
puncture. The wave function for the ID universes, in the true sense of the term (or according 
to traditional definition of the term), is the smooth one point function ipil) = I '^{l) of the loop 
operator. This vanishes in the UV and hence has no arbitrariness in the initial condition. It 
has a growing and then rapidly decaying behavior capturing inflation. 

Let us now argue for the uniqueness of the choice of our trial solution ()6.9|) . As the WdW 
equation ()5.2(i|) is of infinite degree, in principle one needs an infinite dimensional initial condi- 
tion to uniquely choose the solution. Here we argue that a gaussian profile of the initial (/ = 0) 
scale fluctuation 

*(0 + 51) = ^oe""^^'' , (6.20) 
fixes the initial values of all the derivatives 

df--''m = o, ^^^(o) = ^j(rT)^*o. (6.21) 

This provides the infinite dimensional initial condition to uniquely choose the wave function of 
the form \l/(/ < /) ~ e"^'' as the solution of (jSI2EI)- 

Let us now explain the significance of the gaussian profile of the initial scale fluctuation. 
The density perturbation due to the scale fluctuation of the one dimensional universe we are 
considering, is given by 

6p 61 ^{1 + 61) 

j-j = Pi^i-,si = , (6.22) 

where we define Pi^i+si to be the probability of finding an universe of scale / between / and 1 + 61. 
As we know, in an exponentially expanding universe, acts as a cut-off wavelength. The 
vacuum fluctuations with a wavelength greater than freezes as a classical field, producing 
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density perturbations [121 120] • Now using ()6.18p in ()6.22|) 



{l-2CiyX^ 61/1 + 0(51/1)^) , (6.23) 



where \ = H ^. This gives 

SI X' 



I A2 + 2CP ■ 

Computing the quantity ()fi.22|l for / > A = i/^^ we have 



(6.24) 



This imphes a nearly nonflat spectrum of density fluctuation at / > A = H~^{i.e. at big bang). 
For example, for A ~ /^ ^^ (recall / << 1 at the initial moment), in order to achieve this bound 
on the wavelength, 

^ ~ . (6.26) 

This implies a scalar spectral index Ug = 0.96. Thus, presumably one could argue that the 
gaussian initial scale fluctuation profile that uniquely chooses the wave function is not an adhoc 
assumption, rather it physically arises from a special scale dependence of the density perturba- 
tions. 



6.3 The large scale behavior: far future de Sitter 

Let us now consider the large scale limit / > /. We would like to know whether the wave function 
in this region behaves like outgoing plane wave. This would be interesting cosmologically, as 
outgoing (incoming) plane waves in this region would represent expanding (contracting) global 
De-Sitter geometry in the far future [IH]- However, universe at large scale is a classical object. 
So the semiclassical wave function extrapolated to far future is able to capture such expanding 
geometries. The inability of the semiclassical wave function is in capturing the inhomogeneities 
that sets in due to quantum fluctuation in the early universe. 

In the large scale limit I > I, the WdW equation ()5.26|) can be written as 

A cosh ([^) + 5A„,sinh (Z^) - - a] ^(/) = . (6.27) 

Here we have rewritten ()5.2()|1 appropriately for the wave function proportional to a plane wave 
in / space and kept only the contribution from the anomalous dimension that is significant in the 



44 



large scale limit. Clearly the wave function at large scale can assume (outgoing and incoming) 
plane wave solutions of the form \l/(/) = e^*^' with 

dA gdA , , 

Here we have used the fact that around c = 1 fixed point, the betafunctions behaves as (3g ~ 
—g/2, /3m2 ~ 0, with Ng"^ = 1/gst = const. This forces JF to vanish automatically. The 
corresponding spectrum of Dirichlet boundaries at large scale (the one dimensional universes 
with homogenous field) is then given by 

A cosh(i/ A) ± B\n smh{U A) - A = . (6.29) 

This implies the wave function at large scale is an outgoing plane wave (or a far future expanding 
De-Sitter geometry) with the spectrum 

2m7r 2m-7r 



I 9 



-{M^ + n'/W), m,n = 0,±1,±2,... (6.30) 



Again, the non-universal part of the spectrum is only due to Fourier modes of matter and 
factorizes in the wave function. The universal part of the spectrum goes like A^ = '^:S1iML_ 

6.4 The de Sitter minisuperspace? 

Here we will show that suppressing the effect of higher order derivatives with respect to the 
scale in ()5.26|) leads to a de Sitter minisuperspace type of WdW cosmology, namely, a global de 
Sitter in the far past instead of a tunneling wave function for an inflationary universe. The far 
future dynamics of global de Sitter is unchanged, as would be expected. 

For this purpose, let us consider the expansion of the hyperbolic differential operators in 
()5.26|1 using / as the expansion parameter. This is consistent as / << 1 around a c = 1 fixed 
point. We also consider / < / and B/A ~ 0{1/M) 0. The denominator of the fractional 
polynomial in / in the equation ()5.26|) can be expanded in powers of ///. Thus keeping up to 
0{l/lY and second order derivatives in /, (j5.26|) can be rearranged as 

(f T2 ,d 2 Ti 



- + — - ' f 1 - ^) ^1 ^(0 = . (6.31 

dP AP dl n\ A J \ ^ ' ^ 



AP dl p 

This is structurally very similar to a the minisuperspace WdW equation. To make contact with 
all genus minisuperspace WdW equation, one can in principle allow more terms in the potential, 
say term, by truncating the fractional polynomial potential in ()5.26|) at orders higher than 
0(///)^. We identify the factor ordering ambiguity p and the world sheet cosmological constant 
/i to be 
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(6.32) 



with C given by ()6.1Hj) . However, the De-Sitter sign of the cosmological constant 4// < clearly 
shows that the WdW equation derived with all modes of ip turned on {i.e. derived with A^^ 
quantum mechanical degrees of freedom) but truncated to space with quadratic scale derivatives 
is a De-Sitter minisuperspace WdW equation like that of a supercritical theory. This is different 
from the usual AdS minisuperspace WdW equation (|2.4p or (|4.7p of the (non) critical theory 
derived with zero mode truncation of dilaton (f^ [i.e. with quantum mechanical degrees of 
freedom) in a quadratic space of scale derivatives. This indirectly indicates a growth in space- 
time directions. The wave function given by ()6.31|) behaves as the modified Bessel function 

vl/(0~/^(^+^)j_.i(,+i)(2v/^0, (6.33) 

and thus would go like the De-Sitter minisuperspace wave function Jq{\/—^1) at p = —1. This 
represents a global de Sitter in the far past jlHl- 

Recall that the factor ordering ambiguity in the minisuperspace WdW equation (j2.4j) in 
Liouville background is p = —1. This makes the wave function purely Bessel function and 
behaves as plane wave at / — rather than a tunneling wave function. Actually the solutions to 
fl2.4|) are different Bessel or Hankel functions depending on appropriate boundary conditions and 
the sign of /i. They give rise to different classical geometries corresponding to homogenous FRW 
universes. For example, Bessel Kq{2^ — jil) gives rise to AdS geometry (/i > 0) corresponding to 
big bang and big crunch while Bessel Jo(2^— /i/) and Hankel Ho{2y/—fil) can be shown to give 
global de Sitter geometry (yU < 0) in the far past and far future respectively jlHl- These solutions 
are unable to capture any kind of inflationary scenario that introduces large inhomogeneities 
with evolution of time. In our large A^ RG, we see that we have a more general factor ordering 
ambiguity that can admit other values too. For p = —1 we recover the structure of the usual 
minisuperspace WdW in MQM. 

However, the quantum wave function in the / > / regime for Dirichlet boundaries has similar 
behavior as the minisuperspace wave function extrapolated to the large scale. For example 
the I > I solution of ()6.27|1 . outgoing plane waves characterizing a far future de Sitter, can 
be obtained from the large scale extrapolation of the minisuperspace wave function, a Bessel 
function. Thus the far future solutions for universe with homogenous field can be qualitatively 
read from the minisuperspace wave function, presumably this is because the inhomogeneities at 
the large scale are captured more successfully by the dynamics of the nonzero modes of matter 
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that is not dealt by the Dirichlet boundaries. Perhaps one needs to learn this dynamics from 
Neumann boundaries {D-1 branes) in 2D string theory. 

Let us make a brief comment here regarding the magnitude of the de-Sitter cosmological 
constant as captured by large RG. From ()6.32|) we observe 



However, apparently the magnitude is not small due a large H (small g*) around a c = 1 
fixed point, though along the world sheet RG the system flows away from such a fixed point to 
a supercritical regime and we do not know eventually to what value of H (or g) it settles to. 
Nevertheless, as we will see in the last section, nonlocal processes like emission of baby universes 
self-tunes A^s to a small positive value. 

7 The baby universes 

Let us now discuss the nonlinear part of the WdW constraint that connects the different topolog- 
ical sectors in the superspace. The cosmological processes of topology change, such as creation 
of a universe from nothing, or more generally loop splitting due to reconnection of intersecting 
strings giving birth to baby universe [S^ or evolution of cosmic strings connect different 
topological sectors of superspace (each such superspace sector being characterized by all metrics 
having the same topology). This suggests that the Wheeler-de Witt operator could be modified 
by adding a local operator that has matrix elements between different superspace sectors. As 
suggested by Vilenkin this can be schematically written as 



where / is the superspace variable (here the length of the loop) and n labels the topological 
sector (here the number of disconnected loops or the occupation number of the closed strings). 
Considering topology change to be a local process, it can be assumed that A.„„/ (/,/') vanish 
unless / and /' are obtained from each other by changing topological relations at a single point, 
i.e. A„„/ (/,/') 7^ for n' = n ± 1. Knowing the functional form of A„„/ (/,/') needs to deal 
with string interaction vertices at the full quantum level. Even a qualitative knowledge of its 
functional form in the 2D case will be immensely helpful in understanding the topology changing 
amplitudes in the higher dimensions. 




(6.34) 




(7.1) 
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7.1 Nonlinear part of modified WdW equation 

Here we will derive the emergence of the topology changing amplitudes from the nonlinearity 
in the WdW equation captured by large RG . It is argued in [SSI that such a nonlinearity in 
the WdW equation is essential in the matrix model description of 2D gravity. However, there 
seems to be no control over the numbers and sizes of the baby universes emitted or absorbed. 
Here we see that the world sheet RG is very special in this respect that it gives rise to baby 
universes of vanishing size with amplitudes suppressed as the higher order corrections in 0{1/N) 
to WdW equation. This implies that such processes are suppressed at large scale and at weak 
string coupling. The suppression of the amplitudes comes as a virtue of the vanishing size of the 
baby universes. It will be extremely interesting to understand the underlying phenomena that 
enables the world sheet RG in the large language to control the size of the baby universes. 

Let us recall ()5.18|) and consider the multi-trace or wormhole terms represented by the 
ellipses in the equation. In matrix quantum mechanics they are redundant operators. Their 
effects are lower by relative orders of 0{1/N) compared to the leading contribution. To capture 
their corrections to the single topology part of the WdW equation, let us rewrite the expansion 
from as 



S-^Tr- 



1^ 



Tr- 



J.<l>{x) 



d_ 
dV 



g3 



d 



dl"^ 



+. 



Tr- 



(7.2) 



Here Fgi-s are hyperbolic functions given by ()3.17p . Using the expansion in ()7.2|) into ()5.15|) and 
recalling ()5.8|1 . the inclusion of the topology changing amplitudes in the beyond minisuperspace 
WdW can be arranged as 



n^i) + J2 n / 



(7.3) 

where Ti represents the WdW operator acting on a given topological sector. Let us call this 
part single topology WdW operator. Here TC acts on wave function 'i^i(/) corresponding to a 
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single loop and is given by ()5.19|) 

« = - tI) ("m - "'-ii. + ') + ■^-^'^ (4) ^ ('I) ■ <^-^' 

^/'n-s represent n-point functions or the n-loop amplitudes 

^n(/i,/2,../n;ti,t2,..tn) = (Xr ^^lY^^ . . .Tr ) . (7.5) 

To have a flavor about the nature of the matrix elements between the different topological 
sectors, let us write down first few matrix elements as differential operators 

A,,,{1- n = N6{n J2 ^.-^) (^cosh ([-) + i?A„sinh (/-)) , 



m=l 



m=l k=l 



X ( A cosh ( /-^^ + BXr, sinh (l ^ 



dlJ \ dL 

K,,,{l-l\l\t') = -NH{l')5{l")g^ F,2^^(Acosh (/|) +i?A„sinh ([|)) . (7.6) 

They represent branching of a loop into n loops out of which (n — 1) have vanishing radius. 
Thus the {n — 1) loops are emitted (absorbed) as 'baby universes' from (by) the original one. As 
differential operators the dependence on the length of the mother loop and the daughter loops 
factorizes. However, their interacting nature is unveiled by plugging in the scale dependence 
of the n-loop wave function. It will be certainly interesting to solve (17. 3|) by some numeric 
ansatz for the scale dependence of ifjn by and get the scale dependence of the topology changing 
amplitudes. 

To get a clear idea about the relative order of these multi-trace (wormhole) contributions in 
0{1/N), one can absorb the extra dependence in each amplitude by rescaling the length of 
the vanishing loops as V ^ V = Nl' 0. This arranges the multi-trace or wormhole terms in 
relative strength of 0(1/A^) compared to the single trace terms. For example. The correction 
due to Ai^2 begins from a relative 0(1/A^) compared to the single topology part of the WdW 
equation. Similarly, the correction due to A13 is of relative 0(1/A^^) and the correction due 
to Ai^3 begins from a relative 0(1/A^^) compared to the single topology WdW equation. Note 
that, at a given order the contribution from kinetic term of the matrix quantum mechanics has 
much larger contribution than that coming from the quadratic term or the cubic interaction 
vertices. 

It will be extremely interesting to understand the physics in the world sheet RG that controls 
the baby universes to have vanishing size and hence being suppressed by 0(1/A^). In the large 
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RG, from the point of view of the matrix quantum mechanics, this is very closely attached 
to the matrix nature of and hence to the involvement of all the A^^ quantum mechanical 
degrees of freedom in the dynamics. This indicates that perhaps a proper understanding of the 
physics of these A^^ quantum mechanical degrees of freedom in the strong coupling and at small 
scale is a key step towards real cosmology. 



7.2 Self-tuning universe 

Here we will show how the sub-leading contribution to the linear part of the modified WdW 
equation, namely the contribution coming from the formation of one baby universe, self-tunes 
the cosmological constant to a small positive value. As pointed out in [B^j, the mechanism of 
self tuning of the cosmological constant by uncontrollable emission of baby universes [33] has 
been studied extensively in the literature. However, the main problem is to understand the 
issues like the locality of the coupling of the baby universes to the parent universe. As we have 
discussed above, in large A^ world sheet RG a controlled emission of baby universes of vanishing 
size automatically emerges from the nonlinear part of the modified WdW equation ()7.H|1 . These 
events have local coupling with the parent universe by the virtue of the vanishing size of the 
baby universes. Below we will derive the leading contribution of such event from the nonlinear 
part of the modified WdW equation ()7.3|) and its role in self tuning the de Sitter cosmological 
constant ()6.34|1 . 

Considering the leading contribution from the emission of one baby universe, the nonlinear 
modified WdW equation takes the form 



^ ^coth(7rMi?) f d_ 
N 2M J ^'dl 

However, integrating by parts the 0(1/A^) term can be simplified as 



"^^^^'^^^iV ^ 2M / "^^'^^^'^dU ' ^^V2(/,r;a;,t')) =0. (7.7) 



J dl'8{l')^ j dt'MlJ';x,t') = - J dl'e{l') J (7.8) 

where the boundary term 6{l') J dt'ijj2{l,l']x,t') vanishes ?/'2 being a continuous function of /' 
and 0(/') represents the Heaviside step function. Performing the /' integration for /' > (j7.8p 
becomes proportional to 

d f /■* df 

dl'^ing^J dt'Ml,l';x,t') ~ -e-y ^^1, (7.9) 
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where we have assumed all the eigenvalues of (f)(t) are same, given by \(t), and the factor 
of order e°° comes from the /' — oo end of the definite integral. Considering the periodic 
eigenvalue to be cos(^t), such a contribution will modify the de Sitter cosmological constant 
fl6.34|) by modifying the constant A ^ A{1 — a) to be 

A'.. = 2H^(l- = 2H^(l + , (7.10) 



A{l-a)J V 1-cr 

where a is given by 



a = l^coth(7rMi?) r-^e°° 
N2nM ^ 'J \{t') 

^ coth(7rMi?)sec(^t)tan(^t) e°° > 1. (7.11) 



N2nM ' ' ' 
Thus the de Sitter cosmological constant becomes 

A,s^Ks-2H^{l-'-^^^-±^), (7.12) 

2Ce-^ + — 

where it is possible to get — — < 1 for small C as the numerator diverges in inverse power 

in C whereas the denominator diverges exponentially. Thus irrespective of the magnitude of 
H, for small C the sub-leading effect of the emission of a single baby universe can achieve the 
desired self tuning A^^ — ^ A^^ > 0. It will be extremely interesting to capture the phenomena 
from some kind of nonlocal world sheet action (see the world sheet methods in fl^ I26j). 
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